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ABSTRACT 

The exact response of thin shells loaded dynamically is 
extremely complicated and often cannot be expressed in closed 
analytical form. Approximate energy dissipation solution 
techniques are developed herein for rigid-perfectly plastic 
thin shells of arbitrary shape subjected to short duration 
dynamic loads. The effect of finite deflections is considered. 
Approximate results are compared with some exact analytic 
results for simple problems and with some experimental results 
for more complicated problems. 
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CHAPTER 1 



INTRODUCTION 

In recent years realization has grown in the engineering 
comnmnity that structural designs for certain uses can be 
improved by allowing some plastic deformation. 

Elastic-plastic analysis of structures can be quite 
complex. It has been found, however, that elastic deformation 
can reasonably be neglected if plastic deformation plays the 
predominant role in a problem. Analysis has been simplified 
by considering materials to be rigid-plastic. 

Even with use of the rigid-plastic model of material 
behavior, there are few plasticity problems for which exact 
solutions have been found. Great emphasis has been placed 
on finite difference and finite element methods of solution. 

This thesis investigates the dynamic plastic behavior 
of thin, rigid-per f ectly plastic shells and indicates 
approximate solution methods which can be used alone or as 
an early design aid in conjunction with more complete 
numerical methods. 
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CHAPTER 2 



SOLUTION TECHNIQUES FOR RIGID-PLASTIC STRUCTURES 

The stress level at which plastic flow (yielding) begins 
in a material subjected to simple tension is called the yield 
stress or the yield strength, o^. A perfectly plastic 
material can suffer indefinitely large strain once the yield 
stress is reached. A strain hardening material requires 
increased stress if plastic flow is to continue. The yield 
stress of a strain rate sensitive material depends upon the 
rate at which the material is loaded. The material considered 
in this thesis is rigid-perfectly plastic. 

If a solid is subjected to a multiaxial stress field, 
the onset of plastic flo^^; depends on the nine (six 
independent) stress tensor components The relationship 

of the stress tensor components and the yield stress in simple 
tension can be expressed as a yield function 



F(a. .) = 0 
ID 



( 2 . 1 ) 



If plotted in six-dimensional stress space, equation (2.1) 
is a hyper-surface and is called the yield surface. 

Drucker has shown that the yield surface is convex every- 
where and that its shape governs the stress-strain rate 
relationship of the material [2-1] . 
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If the axes are also used as the corresponding 

(plastic strain rate) axes, the strain rate vector must 
be normal to the yield surface at the point corresponding to 
the state of stress. If the stress point is a corner or 
vertex point of the convex yield surface, then the associated 
plastic strain rate vector must be a linear combination of 
the normal vectors to the surface at the singular point. 

The strain rate vector is, therefore, proportional to 
the gradient of the yield function and 



8F(a. .) 

F = A ^3 

^i j 8a . . 



( 2 . 2 ) 



where X is a positive parameter. 

The principle of normality which led to equation (2.2) 
is often called the associated flow rule. 

Many yield functions or yield criteria have been proposed. 
Two which are commonly used and which are in good agreement 
with experimental results are the VonMises yield criterion 
and the Tresca yield criterion [2-2] . 

The VonMises yield criterion assumes that yielding 
under multiaxial stress conditions begins when the distortion 
energy is the same as that at yielding under simple tension. 

The VonMises yield function in terms of principal stresses is 



^[(a^ - + (02 - a^)2 + (a^ - a^)^] - aj 



0 (2.3) 
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The Tresca yield criterion assumes that yielding begins 
when the maximum shear stress is the same as that at yield 
under simple tension. The Tresca yield function can be 
expressed as 

Max{|a^ - a^l/ \o^ - | / | ) - = 0 (2.4) 

If these yield functions are plotted in Haigh-Westergaard 
stress space (principal stress coordinate axes) , the Tresca 
yield surface is found to be a regular hexagonal prism 
inscribed within the VonMises yield surface, a right circular 
cylinder [ 2-2 ] . 

The state of stress at a point in a plate or shell very 
closely approximates a biaxial state of stress (a^ = 0) . 

Figure 2.1 shows the form of the VonMises and Tresca yield 
curves in tv/o-dimensional stress space. 

A perfectly plastic material can suffer unlimited 
deformation once the state of stress reaches any point on 
the yield surface. However, a structure made of such a 
material will not suffer unlimited deformation until enough 
parts of the structure have reached yield to produce a 
collapse mechanism. 

Determination of the exact static load level which will 
cause collapse of a structure (i.e.,the load at which large 
plastic deformation just begins) is often very difficult or 
practically impossible. The static collapse load can be 
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Fig. 2.1 
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bounded, however, by use of the Upper Bound and the Lower 
Bound Theorems of Limit Analysis [2-3, 2-4, 2-5, 2-6]. 

The upper bound theorem states that a structure will 
collapse under the applied load if there is any kinematically 
admissible pattern of plastic deformation (i.e., one which 
satisfies boundary conditions and yields admissible strain 
rates) for which the rate at which the external load does 
work equals or exceeds the internal energy dissipation rate. 

The lower bound theorem states that a structure will not 
collapse (or will just be at the state of collapse) if an 
equilibrium set of stresses can be found which balance the 
applied load and which are everywhere on or within the yield 
surface . 

Both theorems assiame that at the beginning of collapse 
there is no appreciable change in geometry of the structure. 
Results are the same for rigid-plastic material as for 
elastic-plastic material [2-7] . 

The methods of limit analysis have been used successfully 
to find or bound the collapse loads of beams and framed 
structures (2-8, etc.], circular plates [2-9, 2-10, etc.], 
rectangular plates [2-11], shells of revolution [2-12, 2-13, 
2-14, etc.], and non-symmetric shells [2-15]. 

In the case of plates and shells, generalized stresses 
and generalized strains (defined in Chapter 3) are used in 
place of the usual stresses and strains defined at a point. 

The yield function and the associated flow rule are then ex- 
pressed in terms of the generalized stresses and generalized 
strain rates. 
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If the load applied to a structure is in excess of the 
static collapse load, acceleration will result. Actual 
collapse of the structure may be averted if the duration of 
loading is very short. 

Exact solution of such dynamic problems is more difficult 
than exact solution of static collapse problems. The time 
history of the loading process may affect the results, and the 
shape of the deformation may vary during the deformation 
process . 

Few exact solutions have been found for dynamic plasticity 
problems of plates and shells. Some solutions are exact only 
for simplified yield criteria. 

Exact dynamic solutions have been obtained, for example, 
for circular plates [2-16], cylindrical shells [2-17, 2-18], 
and spherical shells [2-19] , under certain conditions of 

9 

loading . 

Stable geometry-change effects caused by finite 
deflections have also been incorporated for some problems 
with beams [2-20, 2-21], plates [2-22, 2-23], and shells of 
revolution [2-24, 2-25]. 

The complexity of dynamic problems in plasticity has led 
to the search for relatively simple approximate methods of 
solution. Techniques have been found which, in principle, 
provide upper and lower displacement bounds for a structure 
which is loaded impulsively [2-26, 2-27]. Unfortunately, 
the upper displacement bound technique requires a lower 
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bound (equilibrium) point load collapse solution, which is 
extremely difficult to find for complex structures. 

An approximate method for describing the dynamic 
behavior of rigid-plastic beams and flat plates was recently 
developed by Jones [2-28] . The method assumes a time- 
independent deformation mode shape, ignoring the motion of 
travelling hinges, and accounts for the effects of finite 
transverse deflections. Comparison with experimental results 
in [2-28] and in [2-23] showed good corelation. 

The major development of this thesis is , an extension of 
Jones' approximate method for beams and plates to the 
dynamic behavior of thin, rigid-plastic shells. 
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CHAPTER 3 



BASIC EQUATIONS FOR AN ARBITRARY THIN SHELL 
Section 1. Assumptions 

The development of certain sets of equations basic to 
the deformation process of an arbitrary shaped shell is out 
lined in this chapter. The equations developed are those 
which relate generalized strain rates to the displacement 
rates of the shell (strain/curvature rate equations) and 
those which relate the generalized stresses and their deri- 
vatives to the external load and the D'Alembert forces (the 
equilibrium equations or equations of motion) . 

Several simplifying assumptions have been made in the 
development. The results are applicable to the wide class 
of problems for which the transverse deformation (perpendi- 
cular to the mid surface of a shell) is the predominant 
deformation. 

The major assumptions are listed below. 

1. The shell is thin; i.e., the smallest radius to 
wall thickness ratio is greater than or equal to 
about ten. 

2. Material is homogeneous and isotropic. 

3. Straight lines perpendicular to the undeformed 
midsurface of the shell remain straight and 
perpendicular to the deformed midsurface. This 
is usually known as the Love-Kirchof f assumption. 
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Displacements in the transverse direction and 
their derivatives are much larger than the 
displacements and derivatives of displacements 
tangential to the shell midsurface. 

5. Transverse displacements may be moderately large, 
but the strains which they produce are infinitesimal 

6. The effects of rotary inertia are small compared to 
those of linear or translational inertia. 

7. The normal stress in the transverse direction is 
small and is assumed to be zero. 

Section 2. Differential Geometry 

In order to find the strain rate, curvature rate, and 
equilibrium equations one must be able to describe the 
geometry of the shell. Therefore, some space will be devoted 
to the differential geometry of a curved surface (the mid- 
surface of a shell, for example). 

The raidsurface of a shell is a two-dimensional space. 
Some information about the surface is intrinsic, that is, it 
can be expressed in terms of the two-dimensional surface 
alone. For a complete representation of the surface, however 
other information is required which relates the surface to 
the three-dimensional space. [3-1] 

Consider the surface element shown in Figure (3.1). 

Each point in the surface can be represented by a three- 
dimensional position vector r (with components r^) with 
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Fig. 3.1 
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respect to the rectangular cartesian coordinates X^. Each 



point can also be represented by its coordinates with respect 
to a curvilinear coordinate system which lies in the surface. 



of mutually orthogonal surface curves which correspond to the 
curves of maximum and minimum normal curvature. The spatial 



It should be noted that Roman indices (i, j, k, etc.) 
will be used to denote spatial tensor variance and will have 
the. range 1 to 3. Surface variance will be denoted by Greek 
indices (a,B,y, etc.) which will have the range 1 to 2. 

Consider two adjacent points, and ^ 2 ’ lying in the 
surface defined by As the two points approach each other 

closely, the length of the spatial differential vector dr 
approaches the differential curve element length ds . 

Thus in the limit as ds approaches zero. 



• ct 

The most simple surface coordinate system, n , is the family 



i i ct 

coordinates can be expressed as X = X (rj ) . 



(ds) ^ = dr • dr 



(3.1) 



By use of the Chain Rule of differentiation 




(3.2) 



a 



and 



dr * dr = r 




(3.3) 



Therefore 



(ds)2 = r^, r^, dn“dn^ 

a 3 



(3.4a) 
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or 

(ds)2 = a dn^'dn^ (3.4b) 

ap 

where 

a ^ = surface metric tensor 

aB. a 8 

• cx 

Since the n are orthogonal, a ^2 ~ ^21 ~ Equation 

(3.4b) is called the first fundamental form of the surface. 

The surface coordinate system also has a set of base 
vectors, a^. These base vectors are tangent to the surface 
in the direction of the corresponding surface coordinate n . 

The surface base vectors are defined as 



a 



a 



3r 

a 

3n 



(3.5) 



The surface base vectors are needed to relate the surface 
to three-dimensional space. This is done by forming the unit 
normal vector to the surface. Since the a^ are tangent to 
the surface, their cross product will be normal to the surface 
and 

a X a 
a 



(3.6) 



with the order of multiplication chosen so that, by convention, 
the direction of n is from the concave side toward the convex 
side . 
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The additional relationship needed to represent a surface 
is that which relates the change in the normal vector, n, to 
change of the position vector, r, as one progresses along 
the surface. This relationship is called the second funda- 
mental form of the surface and is expressed as 



In order to find the elements b one must differentiate 

Olp 

the unit normal vector n with respect to the curvilinear surface 



dn • dr = b , dri'^dri^ 
ab 



( 3 . 7 ) 



where 




coordinates 



Recall that a = — — 
“ 3n“ 



and that r is the position vector 



r = b .r ^ = b . 
1 J 



Since the spatial base vectors, b ^ , for a rectangular cartesian 
system are constants. 




( 3 . 8 ) 



where 



denotes 

a 
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which is the rate of change of the spatial coordinate with 
respect to that of the surface coordinate 

The numerator of equation (3.6) using the orientation of 
Figure (3.1) is 



ai X a 2 = (bj X^) x (b^, X^) 



(3.9a) 



which is expressible in tensor notation as 



a^ X a- = e. X^ b. 
1 2 ijk 1 2 1 



(3.9b) 



where is the alternator symbol. 



e. = 1 
ijk 



for ijk an even permutation of 123 
for ijk an odd permutation of 123 
for two or more indices being equal. 

The denominator of equation (3.6) is the magnitude of the 
numerator. 



e. = -1 
ijk 

e. =0 
X3k 



|a^ X a 2 I = I 1 Ia 2 | sin9 (3.10) 

. ct 

Since the h are specified to be mutually orthogonal, 
sin0 = 1. 

Furthermore, since the magnitude of a vector is the square 
root of the scalar product of the vector with itself. 
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1 



a I = [r^, r^, = /a 

a' a a aa 



(No summation on a) 



Equation (3.10) therefore reduces to 



ai X a^l = 



(3.11) 



and equation (3.6) takes the form 



£. . - b. 

n = - ^3^ A 

^^11 ^ 22 ^^ 



(3.12) 



To find the absolute rate of change of the normal vector 

with respect to the change of position on the curved surface 

one must subtract the apparent change due to the curvature of 
ct 

the n coordinate system. Thus the tensor (or covariant) 
derivatives of n must be used. 



n/a = b^ n^/a (3.13) 

where 

i/ - i .pi jvk 
n /a = n , + F n-* X 

a ]k a 

and Fjj^ is the Christoffel symbol of the second kind for the 
three-dimensional space X^. 

The X^ being used form a rectangular cartesian coordinate 
system and 
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Therefore, the covariant derivative of a space vector 
with respect to the surface coordinates is identically equal 
to the partial derivative of the space vector with respect to 
the surface coordinates. 

That is 

n^/a = n^,^ (3.14) 

Equation (3.7) can be expanded using (3.8) as 

dn • dr = (n/a dn*^) * (b . dn^) (3.15a) 

3 P 

which in turn, by using equation (3.14), may be cast in the 
form 



dn • dr 



n / X Q dri dr] 
a p 



(3.15b) 



Comparing equations (3.15b) and (3.7) shows that b^g may 
be found from 



b 



a3 





(3.16) 



The normal curvature of a surface curve is related to 
the tensors a^^^ and b^^ [3-1] . Passing through every point 
on the surface there are two mutually orthogonal curve segments 
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, whose normal curvatures have extreme values. If the coordinate 



curves 



a 



n 



are aligned with these lines of principal curvature. 



b 



12 




0 



and the principal curvatures are 




and 




(3.17) 



Where R is the radius of curvature of the coordinate 
a 

a 

curve r) . 

In order to fully describe a shell one must bring in a 

third local coordinate normal to the midsurface. This 

direction, Z or ri^, is used to express shell thickness. 

Therefore we shall now consider a family of surfaces 

which are parallel to the shell midsurface. That is, surfaces 
each of which is a constant distance Z away from the midsurface. 

The sense of Z agrees with that of the unit normal vector n. 

The surface metric tensor of a parallel surface is 

designated a q. Since the surface is parallel to the mid- 

CL p 

surface, it has the same principal directions and 



^12 = ^21 = “• 



The nonzero values of a „ are found from the application 

CL p 

of equation (3.1) to the parallel surface at distance Z from 
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the midsurface. The position vector f to a point on the 



parallel surface is 




i . r, i 
r = r + Zn 

and 


(3.18) 




(3.19) 



Proceeding as in equations (3.3) and (3.4) one defines 





(3.20a) 



which, when the indicated multiplication has been carried out, 
becomes 



^ 1 1 ..,1 1 ,,^1 
ao=r, r,-+Zr, n,„+Zr,_ 

aB a B a B B 


4 . 

n m • 2 n / n / « 

'a a 6 




(3.20b) 


Since the position vector component r^ is 


the coordinate X^, 


the second and third terms of equation (3. 
through the use of equation (3.16) as 


,20b) may be expressed 


r, k k 

Zr , n = Zb~ 

'Y o oY 


(3.21) 



29 



The Weingarten formula [3-1] for a unit normal vector 



directed outward states that 



n , = a ' X 

a 6a Y 



(3,22a) 



which permits expression of the third term of equation (3.20b) 
as 



1 1 
^ 'a ^ 'B 



a'5'^ b, = a*'^ a^" b, a 

6a Y X3 y oa A3 YV 



(3.22b) 



since 



X 



i 

Y 



X 



y 



a 



yy 



For the principal. 



orthogonal coordinate system n 



a 



9 



and 





for a = 3 



for a 3- 



Therefore 



1 1 
n , n , 

a a 



aa 



aa 



(no sum on a) 



(3,23) 



Applying equations (3.21) and (3.23) to equation (3.20b) 
yields the result 
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(3.24a) 



a 



aa 



a 

aa 



[1 + 2Z 



b 

aa 



+ (^) 2 ] 
a 

aa 



(no sum on a) 



However 




and, therefore. 



(no sum) 



aa 



= a (1 
aa 



Z ^ 2 
Ra^ 



(no sum on a) 



(3.24b) 



With this result a shell of finite thickness can be 
described in terms of the first and second quadratic forms 
(and tensors) of the shell middle surface and directed 
distance normal to the midsurface. 



Section 3. Strain Rate and Curvature Rate 

The geometric representation of a shell developed in the 
preceding section may now be combined with the definition of 
strain. The approximations of Section 1 are then used to 
develop the generalized strains of a thin shell. 

Consider the local normal direction, the Z coordinate. 

It is also the third coordinate in the right-handed, orthogonal 
curvilinear system n\ r\^ , n ^ = Z. The corresponding metric 

tensor component is a^^ = !• 
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As the shell deforms, let the coordinate system also 
deform such that the coordinates after deformation are 
numerically equal to the coordinates before deformation. In 
general, the metric tensor will change to accommodate the 
deformation. Since the normals to the midsurface remain 
normal during deformation, ^31' ^23 ^32 

remain zero. 

Interest has been limited to finite-deflection with 
infinitesimal strains, and for convenience the undeformed 
shape is used as reference. The description of motion is 
therefore Lagrangian. 

The Lagrangian strain tensor in terms of the three- 
dimensional curvilinear coordinate system and in terms 
of the covariant and contravariant vector displacement 
components and is [3-2] : 







I . / . + y^'^^ y, / . ) 
i/3 k/3 



(3.25) 



where ( )/^ denotes covariant differentiation of ( ) with 

respect to n^. 

Since the local three-dimensional coordinate system 
is orthogonal, 






a 



ii 




and 




(no sum on i) 



(3.26a) 



(3.26b) 
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The form of the covariant differentiations referred to 



in equation 3.25 is 



j/i 



" ^ji 



(3 



and 



u/i j . T^j ^ 

' = + r^. p 



(3 



where ( ),. indicates - and the Christoffel symbol 



mn 



•p 1 ^ \ 

mn * mk,ri ^nk,m ^mn,k'^ 



(3 



i i . a 

But a -^ = a . . = 0 for i j and F is different from zero 
1 j mn 

only if k = Then 



•p ^ __ 1 ^ Sj Sj r ^ ✓v *1 

mn ~ ^ ^mS,,n ^nS,,m ^mn,£,'^ 



(3 



(no sum on £.) 



Furthermore, if Z, m, and n are all different, all three 

terms in the bracket of equation (3.28b) are equal to zero 
£ 

and thus F is zero, 
mn 

I 

Therefore F is different from zero if and only if at 
mn 

least two of the indices Z, m, and n are equal. 

The coordinate curve = Z is a straight line and 



a-- = 1. So a_, . = 0. Also since a._ = 0 for i ^ 3,. 
33 33 , j i3 

a . - . = 0 . 

i3,d 



27a) 



27b) 

is 

28a) 



28b) 
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r 




£ 

The three components of F can therefore be written as 

mn 



mn 





mn 




(3.29) 



Let us look now at one term of the strain tensor E. . 

ID 




Before expanding this representative term, let us invoke 
the assumption that transverse deflections w are much larger 
than the in-surface physical deflections u^ and U 2 » There- 
fore, the only terms retained which involve products of 
displacements (or products which indicate derivatives of 
displacements) are those for which w is involved in each 
factor. Because the Z coordinate is a straight line, 
w^ = w^ = w. With this fact and the approximation noted 
above, rewrite equation (3.30) as 



^11 = ^1/1 ^ '"/I '"/I 



(3.31a) 



which becomes 



E 




(3.31b) 
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after performing the covariant differentiation and 
eliminating the mixed product displacement (because of the 
assumption that transverse displacements are much larger than 
in-surface displacements) . 

Expanding this by using equation (3.29) yields 



; 1 1 



11 



= y 



1,1 



‘ 11,1 



g22 



11,2 



" 11, 3 



■w + 



2 



(3.31c) 



The relationship between the strain tensor and the 

physical strain tensor is derived from the expression [3-2] 



2E. . 
il 



[a. . a. .] 

11 33^ 



-1 = (1+e . . ) (1+e . . ) sin y . . - 1 
^ 11 ' jj ' 'id 



(3.32a) 



(no sum) 

where 6^^ is the amount by which the right angle between 
coordinate curves and changed during deformation. 

For infinitesimal strains, << 1 (no sum) and 

the physical strain tensor can be expressed as 



E. . 

€. . = ^ i 

[a. . a. .]^ 
11 33' 



(no sum) 



(3.32b) 



where e . . + e . . = y • • ^or i j . 
13 31 '13 ^ -* 
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The vector displacements y are related to the physical 



displacements by [3-1] 



u . 
1 




v/rrr 

11 



and 



u 



i 




i 



y 



(no sum on i) 




(3.33a) 



(3.33b) 



and, as noted earlier. 



w = = w. 



(3.33c) 



The physical strain tensor component is therefore 



e = Uj)^! ^ ^11^1 ^ ^11,2 ^2 

2 a^^/5^ 2/1^7 



a 



11,3 



w 



2 




1 hi 



The remaining physical strain tensor components are found in 
the same manner. The nine (six independent) physical strain 
tensor components are 
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HI 



{/s 



11 



11 ^1^ ,1 



^11 



^ 11,1 ^1 ^ ^ 11,2 ^2 ^ 



2/a^ ail 



2/1^ ^11 



^11,3 ^ ^ ^^,1^' 
2 a 



11 



2 a 



11 



(3.34a) 



■22 



(/^22 ^ 2 ^ , 2 



‘22 



^ 22,2 ^2 ^ ^ 22,1 ^1 



2/aJJ ^22 



2/1^ ^22 






2 a 



22 



2 a 



22 



(3.34b) 



e 



33 



= 0 



(3.34c) 



^12 ^21 



*22 






/"A" 

/a 



11 



/t 



22 









^‘^^ll ^22 



(3.34d) 



37 




.1 



(3.34e) 



^13 ^31 



^^11 ^ 1^3 ^ 



^ 11,3 



2/E 



11 



2/a 



11 



2 a 



11 



^23 ^32 






2/a 



22 



2/a 



22 



^ 22,3 ^2 



‘22 



(3.34f) 



Further simplifications can be made through use of the 
Love-Kirchof f assumption that straight lines normal to the 
undeformed midsurface remain straight and normal to the 
deformed midsurface. The in-surface components of displace- 
ment u^ can be expressed in terms of the displacement u^q at 
the midsurface, the slope (u^ Z) q of the rotated normal, and 
the directed distance Z of the point along the normal. That 
is 






(3.35) 



Recall that the surface metric tensor at height Z is 
related to the midsurface metric tensor by equation (3.24b). 



a 

aa 



a (1 
aa 



+ 



^) 
R ^ 
a 



(3.24b) 



(no sum on a) 

The rate of change of this tensor element with respect 
to Z is 

^aa,Z = 2a^^(l + |-) |- (no sum on a) (3.36) 

a a 
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The limitation of thin shells stipulates that 



a 



and, therefore, that 



1 + ^= 1 . 

a 



Equations (3.24b) and (3.36) are simplified to 



a = a 
aa aa 



and 



2a 



aa 



aa , 2 R 



(3.37a) 



(3.37b) 



a 



Since the normals to the midsurface remain normal. 



£i3 £23 0. 



By carrying out the indicated operations of equation (3.34e) 
with the approximations of equation (3.37) and by evaluating 
the results at the raidsurface (2=0) , one finds that 






",i ^ !io 



/a 



11 



and 



<“2,Z>0. ' 



R.^ 



/a 



22 



(3.38a) 



(3.38b) 
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Substitution of equations (3.38) into equation (3.35) 



yields 



u 



1 



u 



10 



+ Z( 








11 



and 



u 

“2 = “20 ^ ^<R 



20 



w 



,2 






22 



(3.39a) 



(3.39b) 



Let us now substitute equations (3.39) along with equations 
(3.37) into the expression for equation (3.34a), 

collecting those terms of which Z is a factor. 



^10^,1 ^11,1^10 , ^11,2^20 , 

^11 L 3 “ , , ■*" 

11 2/a77 a,, 2 a, , /a7^ 



11 11 



11 22 



— + — — 1 + 

^1 2a Z E 

1 2a^^ 



^10 , 



11 



w 



,11 



11 



^11,1 ^10 ^ 



2 (a2_Q_) 



R, 



^ 11,1 ^,1 ^ ^ 11,2 ^20 
2 ( a ^ j ^) ^ ^^ 11*^^22 ^2 



^11,2 ^,2 

^^11 ^22 



(3.40) 



This result and the similar results obtained from 
equations (3.34 b and d) will be expressed in terms of the 
commonly used generalized strains the midsurface normal 

strains, Yq the midsurface shear strain, the bending 
changes of shell curvature, and k ^2 the twisting change in 
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J 

1 




shell curvature, which are defined by the following set of 
equations . 



^11 ^10 ^"^l 



(3.41a) 



^22 ^20 ^*^2 



(3.41b) 



^12 ^21 = ^0 + 2 Z < 



12 



(3.41c) 



Before writing the three strain expressions and the three 
curvature expressions, let us define two metric scale factors 



A and B such that 






A(n\n^) = [a^^(n',ri^)]^ 




(3.42a) 


and 






B(ri^,n^) = [a22 (n ) ] ^ 




(3.42b) 


Proceeding as in equation (3.40) 


and using the notation 


of 


equations (3.42), the (midsurface) strain components and 




(change of) curvature components 


defined by equations (3 


.41) 


are found to be 






_ ^10,1 ^,2 ^20 w 


(w ,)^ 

-f r ^ 


(3.43a) 


10 A AB 


2A^ 


u_. _ B T Ut,. 

_ 20,2 ,1 10 w 


(w_2>^ 


(3.43b) 


^20 B ^ AB ^2 


2B-" 


_A ^10 + B ^20 ^ 

^0 B^A,2A^B^,1 


,1 




AB 


(3.43c) 
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(3.43d) 



= 1/1^10 _ 






1 A' R, 



A ,1 AB 



B 



= l(!i20 _ ^ Ll(!^ _ Ll) 

2 ^ ^2 B \2 AB ^ A ^ 



- u, - w T ^ w _ 

Oi^ = — f 1 0 ' ^ \ j. B / 3 0 _ ' ^ \ 

12 B ^ARj^ ~A^\2 A ^BR2 



(3.43e) 



(3.43f) 



It should be noted that retention of finite (moderately 
large) transverse deflections gives rise to the nonlinear 
term ine^^, E 2 q Yq* These nonlinear terms are missing 

if all deflections are infinitesimal. 

The strain rate and curvature rate equations are 
obtained by differentiating each of equations (3.43) with 
respect to time. In the case of infinitesimal deflections, 
the strain and curvature equations are homogeneous with 
respect to time. If finite deflections are accounted for, 
the nonlinear term in each of the midsurface strains produces 
a nonlinear term in the corresponding strain rate equation. 

The strain rate and curvature rate equations are 



• . ^ ^,2 “20 ^ . ",1 ”,1 
*^10 A AB **1 ^ 



n B U • W W 

• ^ ^ 20,2 . ^,1 ^10 . ^ . .2 /^,2 

20 B AB R2 B^^ 



(3.44a) 



(3.44b) 
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/ 



d 





w 



1 

AB 



(3.44c) 



V = — f— ^ / 20 . , /I f 2 , f 1 >2 

’O B^A\2 a'b\1 AB AB 



u, w 



; = 1(110 - Lil) + ^ 2 .;^ _ 111) 

1 A^ A ^ ,1 AB ^ R 2 B ^ 



• = 1(1^20 

2 B^ R_ B \2 AB ^R- 



1 

- — ^) 
A ' 



u, „ w 



u..„ w 



9i^ = ( 10 — '!■) + ?.(_2_0_ _ /2» 

12 B^AR^ A^ ,2 A^BR^ B^ ,1 



(3.44d) 



(3.44e) 



(3.44f) 



The fact that the strain and curvature expressions 
(and, therefore, the rate equations) derived in this chapter 
contain approximations should be reemphasized. 

There have been many attempts to derive strain and cur- 
vature expressions for thin shells using the Love-Kirchof f 
assiamption or something equivalent, and there are almost as 
many different results as there have been attempts. Differences 
creep into the curvature expressions due to the approximations 
made and due to the stage at which the approximations are 
made. The author has not found any disagreement in the 
expressions for (midsurface) strain. 

Table 3.1 gives a comparison between some representative 
theories. The Table shows what terms, if any, must be added 
to the corresponding curvature expression of this thesis to 
bring it into agreement with the theory referenced. Symbols 
and directions have been changed to agree with those of this 



thesis . 
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Budiansky and Sanders [3-8] defined a pseudo-curvature 
in order to produce a symmetric tensor from which curvature 
could be found. The comparison in Table 3.1 was made with 
the curvature, not the pseudo-curvature. 

The various curvature results would lead to differing 
results for the equilibrium equations as well. Only the 
set of equilibrium equations which is consistent with strain 
rate and curvature rate equations of this thesis will be 
presented. The finite deflections considered do not affect 
the curvature rate expressions. 
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[3-6] 
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Section 4. 



Equilibrium Equations 



Let us now investigate the equilibrium of a differential 
element of a thin shell. In order to ensure that the 
equilibrium equation results will have approximations 
consistent with the rate equations (3.44), the equilibrium 
equations are developed from an energy approach rather than 
from a strictly geometric approach [3-10] . 

In the theory of plates and shells the resultant forces 
and moments which act on the midsurface at an element face 
are normally used instead of the stresses themselves. The 
resultants are defined by a force or moment balance and are 
called generalized stresses. 

A typical shell differential element is shown in 
Figure (3.2) and indicates the positive sense of stresses, 
forces, moments, and pressures used in this thesis. 

Taking a force balance in the ri ^ direction on the face 
perpendicular to ri^ yields 



N Bdn^ = Bdn^ dZ 

^ -h/2 ^ ^2 



A moment balance in the n ^ direction on the n ^ face yields 



R + Z 

M^Bdri^ = c, Bdn^ Z — dZ 

^ -h/2 ^ ^2 
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Fig. 3.2 
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And a force balance in the Z direction on the n ' face yields 



QjBdT]^ 



Bdn^ 

-h/2 



R2+Z 



dZ 



Applying the same technique to the remaining forces 

2 

and moments and using the thin shell assumption that :=■ << If 

X\ 

the definitions of the generalized stresses become 



N = a,dZ 

^ -h/2 ^ 


(3.45a) 




(3.45b) 


^12 " ^21 " ^ 


(3.45c) 


= / O^ZdZ 


(3.45d) 


M 2 = / 02ZdZ 


(3. 45e) 


''12 = -''21 = ^ ^12"^2 


(3.45f) 




(3.45g) 


Q 2 = / T23dZ 


(3.45h) 



The integration limits are the same in each definition. 
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Because of the Love-Kirchof f assumption that normals to 
the midsurface remain normal and straight, ~ ^23 ~ 

Therefore there is no transverse shear strain energy, and 
and Q 2 are considered to be reaction forces necessary for 
transverse equilibrium. 

The affect of the shell thickness has been taken into 
account in the definitions of equation (3.45) and the 
resultant forces and moments are assumed to act on the shell 
midsurface . 

Figure (3.3) illustrates the variation of the generalized 
stresses and reactions of a loaded shell element. 

The principal of virtual velocities states that if any 
stress field which is in equilibrium with the applied loads 
acts in conjunction with any arbitrary, continuous virtual 
velocity field which satisfies the actual boundary conditions, 
then the external work rate (due to reaction forces and applied 
load) is equal to the internal energy dissipation rate. 

Symbolically this is represented as 



W = W 

external internal 



(3.46) 



The internal energy dissipation rate is determined by 
integrating the strain energy rate over the volume of the 
body. This result can also be obtained by integrating the 
product of the generalized stresses with the corresponding 
generalized strain rate over the undeformed midsurface area. 
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Fig. 3.3 
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1 



"int = ^ ” 2^2 ”l2l'0 ^ Vl * «2"2 + “^2’=12''»® 

(3.47a) 



The external work consists of contributions from the 
work done by boundary forces and moments, by the D'Alembert 
forces (inertial forces) , and by the load. 

Figure (3.4) indicates the boundary forces, moments 
and associated displacement. The boundary is a curve and 
the boundary work rate is, therefore a line integral. 



W ^ = 0[N „ 
ext aB 

^boundary 



t 0 + Q n + M .((> ^]dC + 
aB a a aB aB 



/P.u. dA - /yii.u.dS (3.47b) 

S ^ ^ S ^ ^ 



First let us evaluate a representative term of 
equation (3.47a) by substituting equations (3.44) for the 
strain and curvature rates . 

The first term of equation (3.47a) becomes 



u 



/NlCiodS 



= / / [N,- 
12 - 1 - 
n n 



10,1 



+ N 



^20^,2 

L AB 



+ 



N^w 






(Adri^Bdri^) 



(3.48) 
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Fig. 3.4 
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The first and last terms on the right hand side of 

equation (3.47) contain surface derivatives of displacement 

rates. Let us integrate those two terms by parts in order 

to produce a line integral at the boundary which can match 
« 

that found in W . .* * 

ext 



^ 1^10 1 
//[ V '- ' 



N,w ,w , 

' ] ABdn^dn2 = 



f [f (N^Bu 



10,1 



N^Bw ,w ^ 

+ - -A- ’ ^ -)dni] dTi2 = 



N Bw 

/{[N^Bu^q] + 



dri2 - 



//[ (NiB) 




N^Bw , • 

^ A 



drij^ dri2 



where the expression within the {} is to be evaluated at the 
shell boundary. The integral involving that expression can 
be expressed as the line integral 



N , Bw ^ w 

'boundary 



The first term of equation (3.47) can therefore be written 
as 

ot 

*For ease in presentation the n symbols will be typed as 
and the integration paths will be omitted from the //. 
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Nj^ABw 





(3.49) 



Next let us consider the components of motion at a 
shell boundary. Figure (3.5) shows a differential profile 
of a shell midsurface at a boundary before and after defor- 
mation and indicates the displacement rate components on the 
edge = constant. 

The angle is the sum of the angles a and b. 




(3.50a) 



The corresponding angle for the edge t ]^ ~ constant is 



u 



20 



w 



2 




B 



(3.50b) 



From Figure (3.5b) it is seen that 



t^^ = u^Q coscf)^ - w sincj)j^ 



(3.51a) 



n^ = -w coscfjj^ - u^Q sin(f)^ 



(3.51b) 
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Fig. 3.5 
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Similarly, 



'22 



= u 



20 



cosc})^ - w sin (})2 



(3.51c) 



= -w cos(J)^ - U 2 q sin4>2 



(3.51d) 



Using small angle approximation and neglecting mixed 
products of the displacements one finds the components 



hi = '^21 = “lO 



W -W 



(3.52a) 



n 



1 



-w 



(3.52b) 






11 





(3.52c) 



plus five other components which can be obtained by 
permutation of the indices. 

The difference between the external and internal energy 
rates must be zero. If the procedure which led to equation 
(3.49) is applied to all terms of equation (3 . 47a) and the 
results subtracted from the external work rate (3.47b) 
formed with the aid of equations (3.52), the final result 
must equal zero. 

Since the displacement rates are arbitrary, the coefficient 
of each displacement rate must be identically zero. The line 
integral terms, with two exceptions, simply show that the 
internal forces and moments at the boundary are equal to the 
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external forces and moments at the boundary. Since 
at the boundary, the line integral terms which relate to 
the internal moments at the boundary are presumed to be valid 
everywhere within the shell 

(3.53a) 
(3.53b) 

The area integral terms lead to the following results. 



Ql = AB - (V^1 - (^12^^2 - ^\2^,2^ 



^2 = AB tv, 2 - tv, 2 - tV®\l - ^12®, 1^ 



AB(P^ - - V ,1 tN ^ 2''),2 V^l + 



(M^B) T M„B ^ (M^_A) , M^^A - 

N_A .. + ^ 4-^ + ^ ^ = 0 (3.53c) 



12“,2 



R, 



R-, 



R-. 



AB(P2 - yu2o) - N^A^2 - t''l2®^l + V)^2 V,1 + 



m2A) ^2 _ V ,2 ^ tM ^ 2 ^),l ^ V^l ^ Q 



R. 



R. 



R. 



R. 



(3.53d) 



AB(P^ - yw) - 



N^AB N AB N Bw N Aw 

_i - _± + f-J: cJl) + (_± cA) + 

R., ^A \i 'b \2 



(M,B) (M A) 

t^l2'^,2\l tNi2'^,l^,2 t a"^^, 1 t b"^^,2 



M- _A - M., _B , M,A _ 

^ 2«12,12 ' <^',2 

MB 



(3.53e) 
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Retention of the nonlinear, finite deflection, terms in 
the strain rate equations leads to the terms involving the 
transverse deflection w in equation (3.53e), which is the 
equilibrium equation in the Z direction. These terms do not 
appear in the case of infinitesimal deflections. 

The infinitesimal deflection equilibrium equations of 
this thesis have been compared with those of several authors 
for shells of revolution and have been found to be in exact 
agreement. 
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CHAPTER 4 



THIN SHELL MATERIAL RELATIONS 
Section 1. Yield Surfaces 

It was noted in Chapter 2 that the yield surface and 
the associated flow rule for rigid-perfectly plastic thin 
shells are usually described in terms of the generalized 
stresses and generalized strain rates defined in Chapter 3. 

One of the major assumptions of shell theory is that the 
transverse normal stress is negligible, and that the state 
of stress at a point on any parallel surface of the shell is 
essentially biaxial. The state of stress can be expressed, 
then, in terms of two principal, in-surface stresses, 
and The yield surface of a biaxial state of stress 

expressed in terms of the principal stresses was seen to be 
a plane, closed, convex curve (Figure 2.1). 

Substitution of the principal stresses into the 
generalized stress definitions (equations 3.45) produces 
a yield criterion in terms of four principal generalized 
stresses, tv^/o membrane forces (N^ and N 2 ) and two bending 
moments (M^ and M^) . The corresponding yield surface is 
described in four-dimensional space, as is the generalized 
strain-rate vector. A generalized stress yield point is 
reached and generalized plastic strain can begin at a shell 
cross section only when the yield condition for stress at a 
point has been reached for every point of the cross section. 
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That is, plastic flow at a cross section is constrained if 
any level of the cross section remains rigid. 

No mention has been made of the transverse shear forces, 

Q. Since the shell theory used here assumes that there is no 
transverse shear strain, it is assumed that the transverse 
shear stresses and the resultant transverse shear forces 
have no effect on yielding. The transverse shear forces are 
treated as reaction forces necessary for equilibrixam. 

If the principal generalized stresses are not known, 
the yield criterion will be expressed in terms of the six 
independent generalized stresses and - The yield 

surface is then described in six-dimensional generalized 
stress space. 

The class of problems consisting of axisymmetric 
loading and response of thin shells of revolution is one 
for which the principal stress directions correspond to the 
coordinate directions. The four-dimensional generalized stress 
yield surface is, therefore, directly applicable. The yield 
surface based upon the Tresca yield condition was developed 
by Onat and Prager [4-1] . The VonMises yield surface in 
four-dimensional generalized stress space was developed by 
Hodge [4-2]. To the author's knowledge neither yield 
condition has as yet been extended to six-dimensional 
generalized stress space. 

The VonMises yield condition for a shell of revolution 
is considerably more complex for analytical work than the 
Tresca condition. Onat and Prager 's Tresca condition is 
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considered in this thesis to be the exact yield condition. 

For problems in six-dimensional stress space, and often 
in four-dimensional stress space problems as well, approximate 
yield surfaces are used. These simplified yield surfaces may 
be necessary in order to make the resulting equations tract- 
able. Although the results obtained from approximate yield 
conditions will, in general, be in error, the bounding theorems 
of limit analysis allow the exact results to be bounded with 
two approximate yield surfaces of the same shape, one of which 

circumscribes the exact surface and one of which inscribes 
the exact surface [4-3] . 

The two-moment limited interaction yield surface proposed 
by Hodge [4-4] and shown in Figure 4.1, for example, circiim- 
scribes the Tresca yield surface, while a similar surface 
with a scale 0.618 times that of the original inscribes the 
Tresca surface. The results which would be obtained through 
use of the Tresca surface must lie between the results obtained 
through use of the two limited interaction surfaces. 

A two-direction limited interaction yield surface 
suggested by the Tresca yield surface of a curved beam [4-5] 
and shown in Figure 4.2 also circumscribes the exact Tresca 
surface. The inscribing surface also has a scale 0.618 times 
that of the circumscribing surface (Appendix A) . 

Closer bounding may be possible if the stress points of a 
given problem lie only on a small portion of the yield surface. 
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FIG. 4.2 
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Section 2. 



Discontinuities 



In the mathematical theory of rigid-perf ectly plastic 
shells, the possibility of discontinuities in displacement, 
displacement rate, and stress fields exists. It should be 
recognized of course that what may be a discontinuity in the 
rigid-perfectly plastic, midsurface representation of a shell 
would in the actual shell be a narrow zone across which the 
"discontinuous" field quantity varies rapidly. 

Mathematical discontinuities and the associated notion 
of a generalized hinge are useful ideas in the model of an 
ideally plastic shell. Not all discontinuities are acceptable, 
hov;ever, and those which violate the assumptions of the shell 
theory should be identified and eliminated. 

This examination extends the work of Hopkins on flat 
plates [4-7] to shells. 

First let us examine displacements and displacement 
rates . 

The in-surface displacements and displacement rates 

(u - and u a=l,2) may be discontinuous. Fracture is not 
aO aO 

considered in this thesis, and these discontinuities are 
interpreted as severe locaL deformation. 

A discontinuity in the transverse displacement or 
displacement rate (w or w) would correspond to severe trans- 
verse strain. The basic Love-Kirchof f assumption requires 

« 

transverse strain to be zero. Therefore, w and w are assumed 
to be continuous . 
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The continuity of a quantity across a line in a shell 
imposes certain requirements on the derivatives of that 
quantity. 

Figure 4.3a shows a simple curve r(t) in the midsurface 
of a shell element. The s and n directions are tangent and 
normal to the curve, respectively. The curve F is shown 
in Figure 4.3b passing through a differential element of 
length ds and width dn centered on the point P (t) . F divides 
the shell element into two regions as shown, r"*" and R . The 
points P', Pj^, P^, P^, P" , P^, ^2 boundary 

of the differential element. 

Let us define the discontinuity of a field quantity G 
as one crosses the curve F at the point P as 

[G]p = g"^ - g" (4.1) 

where 

g'*’ = lim Gp„ 

PP" ^ 0 

and 

G — lim Gp I 
pp' ^ 0 
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Fig . 4.3 
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(c) 



Fig. 4.3 (con't) 
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Let us suppose that [G] = 0 all along the curve F . Then 




Therefore , 



lim 

ds->-0 




ds 



'll’ = 0 



(4.2) 



(4.3) 



That is, there is no discontinuity across F of the rate of 
change of G along F. 

If we look at the position of F at two instants of time, 
t and t+dt, (Figure 4.3c) 

^^^P(t) " ^^^P(t+dt) " ° 

Since the curve F is assumed to be moving, the change in 
the discontinuity as a function of time is a function of both 
position and time, and leads to a substantial or material 
derivative [4-6]. 



lim 

dt^O 






(4.5) 



where V is the instantaneous velocity of F normal to the 
curve. 
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Equation (4.5) presents the following three possibilities; 



If V 



= 0 



^9t^ 






= 0 

is not necessarily zero 



If V 0 

either 

or 



iff) 



[3G] 



[3G] 



0 



IfiG 

V^9t 



] 



(4.6a) 



(4.6b) 

(4.6c) 



Applying equations (4.3 and 4.6) to the continuous 

• * 

quantities w and w one finds (recalling that w = -^) 



iff' - 0 



(4.7a) 



li' = 0 



(4.7b) 



If V = 0 



dn 






(stationary hinge) 



not necessarily zero 



not necessarily zero 



[w ] - 0 



(4.7c) 

(4.7d) 

(4.7e) 
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If V 7^ 0 



(travelling hinge) 



[ 



Sw -j 

3n^ 



0 



and either 



= [w] = 0 



(4.7f) 



(4.7g) 



or 



>Ih> 




(4.7h) 



The requirements for continuity of forces and moments 
across the line F at the point P are obtained by looking 
at the equilibrium of the differential element of Figure (4.3b). 
This element is redrawn in Figure (4.4) with the membrane and 
transverse shear forces indicated. 

The condition of vertical force equilibrium can be ex- 
pressed as 



( ; - 



^2 



*■2 



/ ) Q cosw dn + ( / 

P ' P ' 

2 1 



^2 

/ ) Q cos5 ds 
p II n 

1 



^2 









P” 



-( / + / ) N sin6 ds - ( / + / ) N simo dn 



'■i 



P” 



n 



'’i 



^■2 
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Fig. 4.4 
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(4.8a) 






^2 




yw) dnds = 0 



If we now allow and to approach Pj^ and P^ and P^ 
to approach P^f dn and d(J) approach zero. Then 



I [Q„] ds = 0 (4.8b) 



Finally, allow P and P. to approach P to find that 

[Q„]p = 0 (4.8c) 

n F 



Similar operations with the force and moment equilibrium 
requirements in the s and n directions yield the following 
results . 



= 



ns P n P ns P 



(4.8d) 



The conclusion, therefore, is that only Q , N and M 

s s s 

may be discontinuous across the curve F. 
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CHAPTER 5 



RESPONSE OF ARBITRARY THIN SHELL TO DYNAMIC LOAD 

Experience has shown that bodies which have been plas- 
tically deformed often have regions with little or no 
permanent strain separated by relatively narrow regions which 
have been severely strained. Although the regions of severe 
deformation have a finite width, they are often narrow enough 
to suggest the appearance of a line. These lines, or hinges, 
are normally thought of in conjunction with severe bending, 
which causes a discontinuity in the slope of the body. It 
was noted in Chapter 4, however, that discontinuities in 
other quantities can also occur in a body. The term "hinge" 
in this thesis will be construed to mean any line across 
which an allowable discontinuity in the stress field or the 
displacement field (or their derivatives) exists. 

Although line hinges, that is, hinges of zero width, 
do not actually occur, the representation of such a hinge 
can be mathematically consistent within the framework of 
the thin shell assumptions so far introduced. 

Since hinges appear to be lines, one would expect the 
dissipation of energy at a line hinge to be represented by 
a line integral. Such was the representation of the external 
work rate at a shell boundary used in Chapter 3 (equation 
3.47b). It should therefore be instructive to look for line 
integral representations of shell energy dissipation expressions. 

If there is no dissipation of energy at the shell boundary 
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(a restriction which will not be maintained for long) , the 
external energy dissipation rate for a shell is the work rate 
of the external forces acting on the shell. 




+ (P^ “ yw)w}dS 



(5.1) 



Each expression within parenthesis may be replaced by an 
expression involving the generalized stresses. The resulting 
equation, which will not be reproduced here, is simply 
obtained by use of equilibrium equations (3.53c, d and a). 

Each term of the resulting expression can be integrated 
by parts over any region within which the term is continuous. 
This process leads to a line integral along the curve bounding 
the region and an area integral within the region. 

Let us carry out the integration of a typical term, the 
first term in the u^^ integrand, over a region of the mid- 
surface bounded by the curve C (Figure 5.1). 



evaluated at the limits. For a given curve C the values 
of and ri 2 not independent, and the limits must be 




//B(N2U^q) ^dn^dri2 



(5.2) 



where {} indicates that the bracketed quantity is to be 

Til 



i 





Fig. 5.1 
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expressed in terms of ri 2 in order to continue integration of 
the first term on the right hand side of equation (5.2). 

This integration process is indicated schematically in 
Figure (5.1). The narrow strip indicates the integration 
over rij^ and the arrows indicate the r \2 integration. The TI 2 
integration can be thought of as 

(Path RQS) (Path RPS) 

The sign of the second term is changed if the integration 
path is reversed to SPR. The entire integration is then in 
the counter-clockwise direction around the closed curve C 
and may be represented as a line integral. A counter- 
clockwise integration path is considered positive by 
convention [5-1] . The same reasoning applied to an integral 
of the form 

leads to a clockwise integration path which, by convention, 
is considered negative. 

The final result for the dissipation term being considered 
is 
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The closed curve C would be chosen as a line across 
which a discontinuity in (N 2 BUj^q) exists. There may be 
several curves of discontinuity for each term of the 
dissipation integral. The contribution from each term is 
the sum of its associated area and line integrals from (and 
around) each region. Not all terms would, in general, have 
discontinuities at the same curves. Each term considered 
might have its own particular regions of continuity and 
hinge lines. Therefore each term should be separately 
summed. For simplicity in presenting the results of converting 
the dissipation integral to the sum of line and area integrals, 
only one summation index will be used, merely to serve as a 
reminder that summation of the contribution from the different 
regions is required. If the procedure described above is 
applied to each term of the dissipation integral which can 
have an allowable discontinuity, the final result is 
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If the external energy dissipation rate expressed in 
terms of external load and displacement rates (equation 5.1) 
were equated to the external dissipation rate expressed in 
teirms of the generalized stresses and the displacement rates 
(equation 5.4), and if the generalized stresses were expressed 
in terms of the displacements and the displacement rates 
through use of the material yield surface and the rule of 
normality, then the displacement field could be found as a 
function of time, in principle. The extreme complexity of 
such an endeavor has precluded exact solution of such a 
problem except for very simple cases. 

In the next chapter the general results of this 
chapter will be specialized to a shell of revolutiorv and an 
approximate solution of the shell response to dynamic load 
will be presented. 

It should be noted that the general dissipation equation 
(5.4) can include hinge lines which coincide with the physical 
boundaries of the shell. Therefore, the external work done 
at a shell boundary can also be accounted for, effectively 
removing the restriction imposed before equation (5.1). 
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CHAPTER 6 



AN APPROXIMATE METHOD FOR DETERMINING THE RESPONSE 
OF A THIN SHELL OF REVOLUTION TO AXISYmETRIC DYNAMIC LOAD 

In Chapter 5 the external energy dissipation rate of a 
general, thin shell of uniform thickness was developed in 
terms of area and line integrals involving the generalized 
stresses, the displacement rates, the transverse displace- 
ment (if deflections are finite) , and their derivatives. 

The line integrals represent mathematically the possibility 
of hinges in a rigid-perfectly plastic material. 

In this chapter the general results indicated in 
Chapter 5 are specialized to the case of a thin shell of 
revolution which is subjected to an axisymmetric dynamic 
load. Neither buckling nor nonsymmetric response is 
considered. 

A commonly used coordinate system for shells of 
revolution is indicated in Figure 6.1 along with the positive 
directions of transverse and meridional load and displacement. 
It is assumed that there are no loads or displacements in the 
circumferential direction. 

Table 6.1 relates the symbols used for general shells 
in previous chapters to the symbols commonly used for this 
axisymmetric class of problems. 
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Fig. 6.1 
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TABLE 6.1 



AXISYMI4ETRIC 


SHELL PROBLEM SYMBOLS 


AND THEIR GENERAL 


SHELL PROBLEM COUNTERPARTS 



Axi symmetric 


General 


Z 


Z 


<}> 


^1 


6 


’^2 






^2 


^2 




A 


R2sin(}) 


B 


u 


o 
1 — 1 
0 


w 


w 




^10 


^9 


^20 


^0 


^0 




^1 




^2 


CD 

-e- 


^12 


Q<}> 


^1' ^1' ^1 


^9' Nq; Qq 


M^; N^; Q2 




^1 


^3 


^3 
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since the geometry of this class of shells does not 
vary with 6, the derivatives of all geometric quantities 
with respect to 9 are identically equal to zero. 

Because of the assumed axisymraetry of the load (and the 
shell response) the derivatives of all stresses and dis- 
placements with respect to 9 are also identically equal to 
zero . 

The non-zero rate equations (3.44a, b, d, e) reduce to 



• , w'w' 



= rT R 



Eg = ^ (UCOt(|) + w) 



• • 



1 . . (w'-u)r; 

= fe- * rT^i 



(6.1a) 

(6.1b) 



(6.1c) 



= 



^ COt(}l 
^1^2 



(u-w' ) 



(6. Id) 



where 



( ) ' 



= iil 

■ d(}) 



The equilibriiim equations (3.53) reduce to 

R^R- sincj)(P. - yii) - N.R^ coscji + (N.R., sin<j)) ' + 

L Z (J) W J- <p Z 

(M R sin4>) ' 

— - MgCOs4) = 0 (6.2a) 
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R^R 2 sin(J) (P^ - yw) - N^R 2 sine}) - NqR^^ sinej) + 



N.R_ sin(}) w' 
(j> 2 



(M.R„ sine})) ' 



1 



(MqCOS(|)) ' = 0 (6.2b) 

U 




(6.2c) 



Noting that the differential midsurface area, dS, is 
dS = (2iTR2sin(})) R^d(}) 

and neglecting external work at the shell boundary, the 

external energy dissipation rate can be expressed in terms 
of external loads as 



The dissipation rate in terms of area integrals and 
line integrals may be obtained from the general results 
(equation 5.4) or by applying the procedure of Chapter 5 to 
equation (6.3). Application of the procedure to a portion 
of the dissipation integral is demonstrated in this chapter. 

Symmetry limits hinge curves in the rigid perfectly- 
plastic material to circles ((J) = constant), and the regions 
of continuity for each term of the dissipation integral are 




(6.3) 
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therefore each bounded by two circles. Since there is no 
variation with 0, the problem is essentially one-dimensional. 

Regions in which a typical term of the dissipation 
integrand (expressed in terms of the generalized stresses) 
is continuous are indicated in Figure 6.2. 

Integration by parts of such a typical term gives 



/ (-N R 2 sin(j)) ' u d(j) = S/N , R 2 sin(J)u'd(j) - Z{N R 2 Sin(J)u} (6.4) 
<P n n 



where { } indicates that the expression is to be evaluated 
between the upper and lower limits corresponding to the 
nth region. 

Indicating the value of (N^R 2 sin<j)u) evaluated by 

approaching i> from below as ( ) and from above as ( )"*’, 
m m . m 

the second term on the right hand side of equation (6.4) can 
be expressed as 



- Z {N R Sin4.i} = - Z [( >; - ( );i (6.5) 

n ^ m 

If the summation over m is expanded and the terms 
grouped according to the pLastlc hinge index m, this 
becomes 



- Z {N R sin(|)u} = -[(),-( )^] -[()“-( ) 2 ] - 
n ^ 

= + Z [N , R_sin(j)u] (6.6) 

2 ^ m 

ra ^ 
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Fig. 6.2 
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where indicates the discontinuity in F across the hinge 

circle (p = <j>^. That is, [F] = ” F . 

m m m m 

If attention is restricted to shells with smooth 

generators, geometrical factors such as R 2 sin<}) are not 

discontinuous. Let the notation (G) indicate that the 

m 

function G is to be evaluated at 6 . Then equation (6.4) 

m 

can be written as 

/(-N R_ sine}))' ud(}) = 

<P n ^ 

+ Z(R^ siTuP)^ [N u]^ (6.7) 

m ^ 

Investigation of discontinuities across the hinge 
circles based on the assumption that the w and w are 
continuous everywhere and based on the interpretation that a 
mathematical discontinuity in u merely represents a very 
large local strain yields the following results (c.f. 

Chapter 4) . 

[M^l = [N^] = [Q^] = 0 (6.3a) 

[w'] may be different from zero at a 

stationary hinge. (6.8b) 

[w'l may be different from zero at a 

travelling hinge. (6.8c) 
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that 



Since [M,] and [Q. ] are both zero, equation (6.2c) shows 
(j) (f) 



R- cot({) 



(6.9) 



With the discontinuity results of equations ( 6 . 8 ) and 
(6.9) and the notation of equation (6.7), the external 
dissipation rate expression for the axisymmetric case may 
be written as 



•M 



5 ^= -/,{ZN 'R„ sin(|) + E (^) ' R» sin(J> + I ■ '^ 2 — ]ud(|) 



M.RtR., sin<J) 
4> 1 2 



27T 



(j,^r’e ^'2 



g 1 



m 1 



M,R_ sincj) 

+/^{EN,R_ sincj) - EN-R_ sin(J) + Z-^-- 

(p©2 .62 _ R, 

^ n ^ 1 ml 



M R sin4> . 

}u'd(j) +/^{N^R 2 + NgR^}sin(J) wd<J> 



R'R_ sin<t) 

<” 4 ,“ + »4,> - R, 



R- sin(J) 

- M-COS(i) }w' d(b -/, {E — (N.w + M, ) }w"d(J) 

6 K ^1 



M 



+ E{r sin(J)(N. + 5 ^)} [u]_ - Z(R sin<l)) [NqU] 

q2 <PR3^qqp2 p6p 



R»sin<J) 

y 1 ^ < 1 



R,sin(J){N^w + M^} . 



-) [W] 

K K 
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( 6 . 10 ) 



R_sin(})N . w 



It should be noted that the nonlinear terms in w appear 
only if finite deflections are considered. One of these, the 
final term of equation (6.10), can be shown to be zero. 

Continuity considerations revealed that [w’j could be 

different from zero at a stationary hinge only. Now [w'] ^ 0 

• • 

implies that w" and w" are infinite and, therefore, that 

is infinite. For the exact Tresca yield surface of shells 

of revolution (Onat and Prager [6-1]), normality requires 

• • 

the ratio <./£, to be finite unless N, is zero. Finite 
<j)' (J) 

deflections give rise to membrane forces and, therefore, 
e, would also be infinite. If the shell stretched at a 
hinge with an infinite rate for a finite period of time, 
infinite thinning would occur. The finite deflection term 
involving [w'] is therefore zero. 

Solution of equation (6.10) in terms of displacements 
as functions of time and position is a very complex problem, 
and is made more so because of the presence of travelling 
hinges . 

In order to simplify the problem, let us restrict our 
attention to problems with transverse loading only. 

For the remainder of this chapter it will be assumed that 
= 0 and that ii = 0. Therefore, the integrals and hinge 
terms involving u are identically zero. 
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Another approximation, suggested by Sawczuk [6-2] and by 
Jones [6-3] for the static and dynamic response, respectively, 
of flat plates, is also used. 

Although the discontinuity [w'] can only exist at a 
travelling hinge, let us treat such a hinge as if it is in 
fact stationary and formed at the same location at which a 
static collapse hinge would form. Identical displacement 
profiles for static collapse and for dynamic response have 
been found for cases such as square plates [6-4] , circular 
plates [6-5] , and cylindrical shells [6-6] , if the applied 
load does not greatly exceed the static collapse load. 

If the external dissipation rate expression for the 
transverse load problem in terms of generalized stresses is 
equated to the dissipation rate expression in terms of 
external load, the result is 



E • 

sincj) + sin<j)}wdd) + 

27 r(J)(j )2 91 tt 



R' R„ sin(j) 

Rl (V ■ '”*"2 I7’ 



(}) (j) (j) 2 



r fk M 1 % J ± v’T -R ^ SI rii d) 

W*d(J) - Zf. 2 /x-r . T./T N 



'(j)” 



R 2 sincj) 



Z{— (N,w + M, ) } [w'] = / R-R_ sinct) (P,-yw) wd(}) (6.11) 

^ Rj^ (p (p K < (pi2 '3 



where boundary loading has been omitted for simplicity. 
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To use equation (6.11) in estimating the permanent 
deflection of a shell of revolution subjected to an 
axisymmetric , transverse dynamic load, one begins by 
assuming a kinematically admissible velocity profile 
whose shape is time independent. The static collapse profile 
is suggested, if known. 

w(4),t) = Wq (t) ‘f ((J)) (6.12) 

} 

where f{(})) is the deformation velocity mode and w^ (t) is 
the velocity mode amplitude. 

The velocity profile of the assumed mode of deformation 
is used with equations (6.1) to determine the strain and 
curvature rates. 

The rate equations are used in conjunction with the 
assumed yield surface and its associated flow rule 
(normality) in order to determine the values of the 
corresponding generalized stresses. 

The expressions derived for the generalized stresses 
are entered into equation (6.11) and the indicated 
summations and integrations performed. The result, after 
factoring out common coefficients, is a second order 
differential equation in w^ and time. The equation will be 
nonlinear if finite deflections have been retained. 
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The exact formulation of this differential equation is 
dependent upon several factors such as shell geometry, assumed 
deformation mode, and assumed yield surface. Specific 
examples are given in Chapter 8. 

The method of solving the differential equation depends 
upon its formulation and upon whether or not the pressure 
pulse shape is known. If the shape of the pulse is not 
known or if it is a high pressure, short duration pulse, the 
problem can be treated as impulsive loading. 

If the time history of the pressure pulse is to be used, 
the initial conditions are 

WQ(t=0) = WQ(t=0) = 0. 

If, however, impulsive loading is assumed, the initial 
conditions are 

WQ(t=0) = 0 

w„(t=0) . vVg 

where is the actual impulsive velocity at the location of 
Wq, and V is the Martin and Symonds mode factor for optimum 
matching of the actual impulsive problem with the mode 
approximation (time independent profile) [6-7]. Equation (6.11) 
reduces to the flat plate approximation of Jones [6-3] . 



93 



i 




CHAPTER 7 



APPROXIMATE RESPONSE OF NON-SYMMETRIC SHELL TO DYNAMIC LOAD 

The solution technique of Chapter 5 can be applied, in 
principle, to any shell which satisfies the usual thin shell 
theory assumptions. In this chapter, however, an approximate 
technique is investigated which is not derived from the 
previous results. The technique for dynamic loads developed 
here is a simple extension of a static collapse method pro- 
posed by Janas for the limit analysis of non-symmetric shells 
[7-11]. Finite deflections have not been included. 

The basic assumption used in this approach is that the 
body deforms as several rigid regions separated by narrow 
plastic hinges. The plastic hinges are viewed mathematically 
as curves. The results of severe membrane strain, twisting, 
and bending rates are viewed in terms of discontinuities in 
velocities and relative rotation rates between rigid regions 
at the hinge curves. All internal energy dissipation occurs 
in the plastic hinges, and, according to the principle of 
virtual velocities, the external work rate must equal the 
internal energy dissipation rate. The dissipation rate can 
be expressed as a function of the hinge discontinuities. 

Consider a hinge separating rigid regions I and II to 
have a small width X, as shown in Figure 7.1. Because trans- 
verse shear strain is assumed to be zero, there can be no 
discontinuity in the transverse displacement rate w. The 
only displacement rate discontinuities 
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Fig. 7.1 
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then are ^^20^ correspon«5ing to the shell coordinate 

axes rij^ and 

The velocity discontinuity components normal to and 
tangent to the hinge curve C are seen to be 



[Unl = 




cosB 


'“io> 


sinB 


(7.1a) 


[up = 


l“iol 


cosB 


0 

CN 

1 


sinB 


(7.1b) 



where B is the angle between the normal to the hinge curve 
and the coordinate axis Ti 2 * 

Regions I and II are rigid and the length of the hinge 
curve can not change. Therefore, 



e 



t 



0 



(7.2a) 



If one assumes that the displacement rates vary linearly 
across the width. 



and 






(7.2b) 



(7.2c) 
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By transformation of axes it can be shown that 



£ 



10 



e 



20 



>“ 10 > 



^^20^ 



sin3 



cos 3 



(7.3a) 

(7.3b) 



^0 = 



[u^q] cos 3 + fi^20^ sine> 
Y 



(7.3c) 



By similar reasoning one finds the curvature rates across 
the plastic hinge. 



• 1 

= — sin3 (7.4a) 

• ^2 

<2 = — cos3 (7.4b) 

cos3 + sin3 

2<i 2 = -i j— ^ (7.40 



where is the component of the relative rotation rate 
between rigid regions I and II. 

The internal dissipation rate is 



°i g ^^1 ^10 ^2 ^20 ^12 ^0 + + 

hinge 



”2 ^2 “12 ''l2> 



(7.5) 
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If it is assumed that the generalized stress values are 
uniform across the hinge zone, and if it is noted that 



dS = XdC 



then equation (7.5) can be written in the form 

= / {N^[u^q] sinB + N2[u2q] cos 3 + cos3 + 

C 

^12^^20^ sin3 + sin3 + M 2 ^2 + 

2M^2 cos3 + 2M^2 ^2 (7.6) 

A limited interaction yield surface is assumed which 
reduces to the two moment limited interaction yield surface 
(c.f. Chapter 4) for N ^2 = ^^2 “ It is assumed that 

membrane shear strain can occur only if N ^2 = 7 and that 
twisting can occur only if M ^^2 = 7 M^. The factor of one 
half arises from the fact that the yield stress for shear 
is one half the yield stress in simple tension. 

The limited interaction yield surface is shown in 
Figure 7.2. 

The scale of the similar yield surface which would 
inscribe the exact six-dimensional Tresca yield surface is not 
known because the six-dimensional Tresca yield surface has not 
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FIG. 7.2 
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been developed. The scale of the inscribing surface can be 
no larger than 0.618, which is the scale of the two moment 
limited interaction surface inscribed in the Tresca surface 
in four-dimensional stress space. The value 0.61S is used in 
this thesis as the best available scale factor. 

The solution to a particular problem would be carried 
out as follows. 

First a collapse mode is postulated which involves rigid 
body motion of one or^ more regions. 

The rigid body equations of motion for each region are 
examined at each hinge curve to determine the velocity dis- 
continuities and the relative rotation rates between regions. 
The condition of continuity of transverse velocity w may be 
used to determine the equation of one or more hinge curves. 

All velocities may be expressed in terms of a characteristic 
transverse velocity w^. 

The strain rates and curvature rates are used to determine 
the corresponding generalized stresses. 

The internal dissipation rate (equation 7.6) is equated 
to the external dissipation rate (equation 5.1). The resultant 
second order differential equation in Wq and time can then be 
solved for the permanent deformation w^. 

The procedure indicated in this chapter has been used to 
solve the problem of a cylindrical panel loaded impulsively 
over a portion of the shell inner surface. The ror.ults are 
reported in the next chapter. 
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Extension of this technique to include finite deflections 
is not conceptually difficult. The major difference is that 
the path of a point in the midsurface of a shell region under- 
going rigid body rotation must be recognized as the arc of a 
circle rather than assumed to be the tangent to the arc. 

Rigid body motion cannot be limited, in general, such 
that transverse displacements are much larger than in-surface 
displacements. Non-linear terms involving all of the dis- 
placements could therefore be expected to arise in finite 
deflection problems. 
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CHAPTER 8 



APPLICATION OF GENERAL METHODS TO PARTICULAR CASES 



Section 1; Complete Spherical Shell Subjected to Unifonnlv 



As a first example let us examine the case of a complete 
spherical shell of radius R and thickness h subjected to a 
unifoinnly distributed, exponential internal pressure, p^. 
Comparison will be made with the infinitesimal deflection, 
two-moment limited interaction (Figure 4.1) solution of 
Sankaranarayanan [8.1]. 

The deformation field for this trivial problem is one of 
uniform radial expansion. 

The velocity field is 



Distributed, Exponentially Decaying Pressure 



w((}),t) = Wq (t) - 0 



( 8 . 1 ) 



There is no bending (M^ = = 0); therefore, normality 



requires that 



N 



N 



0 



N 



0 



( 8 . 2 ) 



where 




Dissipation equation (6.11) reduces to 
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2L 

2 



2 / 2Rsini}) di}) = 



0 



!L 

2 

2 / R^sinij) (p 
0 



3 - (8*3) 



Letting 



Pj = Po®*’"''’''’ 



and evaluating equation (8.3), one finds that 



= ^0®'"''"“ -'■o 



(8.4) 



2N 

where P = — 
c R 



= the uniform static collapse pressure of a 
complete spherical shell 



and tQ is the decay time constant. 
Integrating, one finds that 






(8.5) 



since 



WQ(t=0) = 0. 



Motion will cease at time t^, obtained from the positive 
root of 



-P t e ^f'^^0 - P t + P t 
0 0 ^ c f 0 0 



= 0 



( 8 . 6 ) 
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Integrating equation (8.5) with the initial condition 
WQ(t=0) = 0 and evaluating at t=t^ yields the final dis- 
placement 



"Of = ^ 



p_t 



C f , 



(8.7a) 



or through use of equation (8.6), 



"Of = ^<'’o-^c>Vo - 



(8.7b) 



These results are in complete agreement with those of 
reference [8.1J. 



Section 2; Clamped Circular Cylinder Subjected to Uniform 
Impulsive Load 

Consider one bay of a rigid-plastic cylinder which is 
reinforced periodically with rigid circumferential stiffeners. 
Because of symmetry the boundary conditions at each stiffener 
are those of the clamped condition. 

Hodge investigated the response of such a shell section 
for infinitesimal deflections when subjected to axisymmetric 
dynamic pressures [8-2]. 

Jones considered the response of the same cylinder when 
subjected to an axisymmetric impulsive load with and without 
the effect of finite transverse deflections [8-3] . This 
solution included travelling hinges and for the infinitesimal 
case is much simpler than the corresponding dynamic pressure 
pulse analysis of Hodge [8-2]. 
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In reducing the strain rate equations (6.1), equilibrium 
equations (6.2), and the energy dissipation rate equation (6.11) 
to the case of a circular cylinder with transverse velocity 
field, it should be noted that ip is tt/2 radians, is 
infinite, and d(f) is zero. The product Rj^d(J) is, however, 
the differential longitudinal element dx. Using the 
notation 



( )* 



= ( ) ' = < 3 ( ) 

Rj^ dx 



and assuming u to be zero, the strain rate equations (6.1), 
take the form 



£ = w*w* 

X 



e 



0 



w 

R 



X 



-w 



** 



0 



0 



(8.8a) 

(8.8b) 

(8.8c) 



The consistent set of equilibrium equaticms (6.2) reduces 
to 

Q = -M* (8.9a) 

X X 

(p - yu) = N* = 0 (8.9b) 

IK ^ 

^0 

(p^ - yw) = ^ - (N^w*)* - 

when Mq is assumed to be zero. 



* 

X 



(8.9c) 



Equation (8.9b) is a consequence of the assumption 
that u = 0 when p =0. 

The approximate energy dissipation equation (6.11) 
reduces to the form 



D 



E 



2ttR 



Z/{-| N*ww* 

i ^ X 



(N w + M )w**}dx 

X x 



- E (N w + M ) . [w*] . = /(p, - yw)wdx (8.10) 

• X XT I J 

3 J J 

where the summations are taken over i continuous regions and 

over j hinges. 

The velocity profile is assumed to be symmetric about 
the midspan of the cylinder (Figure 8.1). 



w (x, t) = Wq ( t) ^ 



ft < < T 

for 0 - X - L 



( 8 . 11 ) 



where w^ is the velocity at x = L. 

Using the velocity profile (8.11), symmetry about x 
and the equilibrium result that N* = 0, equation (8.10) 
to 



L RL 



w. 



-2w, 



dx 



- ^'Vx=0 -T - 



(N w 

X 



M ) ^ 

X x=L 



^) = 



L MW»x W-X 

2 <P3 - ^ 



= L, 
reduces 



( 8 . 12 ) 
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Fig. 8.1 
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From the two-direction limited interaction yield surface 
(Figure 4.2) it is seen that, since = 0, = N,. for 

D u 0 

w > 0 . 

If one assumes that the central hinge is twice the 
length of the end hinges and that the longitudinal strain is 
concentrated equally in each of the three hinges, it can be 
shown that the ratio of the average longitudinal strain rate 
to the average curvature rate in the hinges is 



e avg. w 

= - 2^ (8.13) 

avg . 

where (+) is used at the clamped ends and (-) at x=L when 
w > 0 . 

The equation of the x-direction portion of the yield 
surface is 



M = ±M- [1 

X 0 



N 

o 



By normality 



(8.14) 




Therefore , 




2N M 




0 



2NxN^h 



(8.15a) 



(8.15b) 
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from which 



N = N- 
X Oh 



(8.15c) 



However, N cannot exceed N. and the final result is 
X u 



w 

N = N- 
X Oh 



w 

. 0 < , 
for r— - 1 
n 



N = N- 
X 0 



f ^0 > , 

for - 1 



M = ±M„[1 
X 0 






f ^ 1 

for jT - 1 



M = 0 

X 



w. 



^ 0 > T 

for 7— - 1 

h 



(8.16a) 

(8.16b) 

(8.16c) 

(8.16d) 



where (+) is used at the central hinge when w > 0. 

These results are the same as those of [8-3] for a 
clamped beam. 

If equations (8.16) are substituted into equation (8.12) 

along with the fact that N_ = N. and the spatial integrations 

b 0 

are then performed, the results, after some simplification, 
are 



2M 

W 



0 



+ 




(1 + 





PWq 

3 



, . , 0 < , 
while r— - 1 
h 



(8.17a) 



and then become 

2Mq 4M^ W(, p3 

Rh h 2 3 



when 



w. 



h 



exceeds 1. 



(8.17b) 
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It should be noted that there are no longitudinal 

membrane stress and no longitudinal membrane strain rate if 

the transverse deflections are considered to be infinitesimal. 

The kinematically admissible velocity profile (8.11) would 

require M to be +M and at the central hinge and at the 

X (J u 

boundary hinges, respectively. 

From equation (8.9c) it can be verified that the 
generalized stress field 



and 






4x 2x ^ 



is statically admissible. 



The upper and lower bound static collapse pressures, P^, 



coincide . 



P 

c 





(8.18a) 



or 



P 

c 





or 



P 

c 




(2 + c^) 



where 




P is the same as equation (11) of [8-2] . 
c 



(8.18b) 



(8.18c) 
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The dissipation equation (8.12) yields the following 



results for dynamic infinitesimal behavior 

“o , _ P3 ““o 

Rh 2 “ 3 



(8.19) 



Equation (8.19) is an ordinary linear second order 

differential equation in Wq ( t) . Equation (8.17) is a pair 

of ordinary second order differential equations, one of which 

w 

is nonlinear and governs the motion until ^ reaches unity, 
and the other of which is linear and governs the motion if 
and when exceeds unity. 

For the case of impulsive loading, p^ is considered to 
be zero. The impulse delivered to the shell is accounted for 
by the initial conditions of the mode approximation. 

WQ(t=0) = 0 (8.20a) 

and 

WQ(t=0) = Wq (8.20b) 

where v is the Martin and Symonds [8-4] mode factor (defined 
below) and is the characteristic initial velocity deter- 
mined from the applied impulse by the conservation of 
momentum. 

[ /pdt = pVq ] 
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Martin and Symonds showed for infinitesimal deflections 
of an impulsively loaded rigid-perfectly plastic structure 
that the best agreement between the physical deformation 
process and a mode approximation (time-independent velocity 
profile) results when the characteristic mode initial velocity 
Wq is equal to a certain multiple of the exact characteristic 
impulsive initial velocity [8-4]. It is assumed herein that 
the same factor is valid for finite deflection problems. The 
factor, V, is shown to be [8-4] 



/ uu. V. dS 

§ ( 8 . 21 ) 

• • 

/ uu. u. dS 
S ^ ^ 



where y is m.ass per unit area, u^ are the mode approximation 
initial velocity profiles, u^ = V^f^ (x) , and are the actual 
impulsive initial velocities = Vq g^ (x) . 

For this particular problem, only u^ and are different 
from zero. 



and 



u., = 



V. - 
0 L 



V.. = V, 



Substituting these values into equation (8.21), the mode 
factor is found to be v = 1.5. 
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The motion equations (8.19) and (8.17) can be rewritten 



in the following form for the impulsive case 



w 



0 



3? 
c 

2v 



for infinitesimal deflections 

(8.22a) 



or, for finite deflections 



6M- 3P 

+ -r^ -I- ^ 

0 yL“ 2]i 



= 0 



while 






and 

« + !o 

0 yL' h 2y yL^ yRh 

when 

T— - exceeds 1 
h 



(8.22b) 



(8.22c) 



where is given by equation (8.18) and the initial conditions 
for (8.22 a and b) are given by equations (8.20). 

The infinitesimal problem is easily solved by integrating 
equation (8.22a) twice and recognizing that motion ceases 
when Wq = 0. 

The ratio of the final central deflection w^^ to the 
shell thickness h is 
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(8.23a) 




h leM^h (Rh+L^) 



or 



h 2 (2+c^) 



(8.23b) 




where X is the loading parameter . 

4M^h 

Solution of the finite deflection problem is more 
involved. Solution of the nonlinear equation (8.22b) with 
its initial conditions (8.20) is sufficient so long as 



If the solution yields 



then equation (8.22c) must be solved for the second stage of 
motion. The initial conditions for the second stage of motion 
are that the displacement equals h and that the velocity 
equals the velocity at the end of the first stage of motion. 

Let us write the first stage finite deflection equation 
in the form 




3P 

2 " 2y 



c 



3Mq (2+c^) 






(8.24) 
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Note by use of the Chain Rule that 



.. _ • '*”0 
^0 ^0 dw^' 



Substitute this expression into equation (8.24) and 
integrate between the displacement limits 0 and to find 



B w ^ 

i {wJCWq) - Wq(0)} = - ( —3 - + ^2 



(8.25) 



Recalling that 



Wq(0) = Wg 



and that Wq is maximum when Wq 



= 0 , we see that 



38„ 

w^ + 3 w 

Om B, 



3vWJ 



Om 2 B. 



= 0 



(8.26) 



Using the definitions of Bj^/ ^2' ^c' value of v 

for this particular problem along with the formula for solution 
of a cubic equation [8-5], it can be shown that 



^ - t|^ <1^ 



^ '1^ 



2 +c 



2 — 7 



) + (|^)" 



(8.27) 



(First stage of motion) 
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calculated from 



Now, 

W- 

Of _ Om 
h ~ h 

if less than unity. 

w 

If, on the other hand, the value of r 

h 

equation (8.27) is greater than unity, find from equation 

""o 

(8.25) that the velocity when = 1 is 

28 h^ 1 

Vi = WQ(h) = Wq(t=0) = [v^Vq - 282 ^]^ (8.28) 



where t is time of second stage motion (i.e., x=0 when WQ=h) . 

For this particular problem equation (8.28) can be 
expressed as 



Wq (h) = 



V7- (t = 0) = 



- 



2+3 (2-i-c^) 
2X 



(8.29) 



The finite deflection, second stage equation (8.22c) 
can be written in the form 






"0 



64 - 0 



(8.30) 



where 



B 



3 



12M 



0 



yL^h 



and 
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6M, 



6 . = 



4 yRh 



with initial conditions 



Wq (t= 0) = h 



and 



Wq(t= 0) is given by (8.29). 



The solution is 



V, 



6 



6 



~ T~ sin fx + (h+^— ) cos fx - 7 ^ 



6 . 



B. 



(8 



and the velocity is 



Wq (x ) ~ cos f 



^4 

- f (h+^; — ) sin f X 
^3 



(8 



where 



f = /Bj 



In terms of this particular problem, the second stage 
solution can be expressed as 



V:', /3> . (8+3cM . ^ iici 



cj 

4 



.31a) 



.31b) 



h 
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cos f X 



(8.32a) 



and 



Wq(t) 




( 8+3c^ ) 
6 



cos ft 



(4+cMf 

4 



sin 



f T 



(8.32b) 



If the shell does enter the second stage of motion, the 
final displacement-thickness ratio is 





(8+3cM 

6 



4+c ^ 

sin fi^ + cos fx^ 




(8.33) 



where x^ is the first positive root obtained from setting the 
second stage velocity equal to zero. 

That is. 



fx 



1 



arctan 



/3X (8+3c^) 

/ 4 " 6 



4+c^ 



(8.34) 



For the simply supported case with finite deflections it 
may be shown that 



N 



2N — 
0 h 



and 



N 



X 



N 



0 



for 



^ < 1 
h 2 



. > 1 
for ^ 



(8.35a) 



(8.35b) 
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The factor in equation (8.24) becomes 



12M^ 

= -5- 



and $2 becomes 






The 



3Mq (1+c^) 
yL^ 

final displacement-thickness ratio becomes 



for 



or 



h ^16 



+ [ 



16 



w 

h 



0 < 1 
2 



+ 




(— ) ^ 




(^)2 



] 



1 

T 



^0 ^ /TX (4+3c^) 
h / 4 12 



sin 




+ 



2+c^ 

4 



cos fx^ 



c 

4 



for 



w 

h 



0 > 1 
2 



4 /3X 4+3c^ 

where fx^ = arctan { — } 

^ 2+c^ 



(8.36a) 



(8.36b) 



(8.37a) 



(8.37b) 
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Equation (8.23) 

Jones (infinitesimal) [8-3] 

Hodge (rectangular pulse lOX static 
collapse pressure) [8-3] 

Equations (8.27 & 8.33) 

Jones (finite deflection) [8-3] 



Fig . 8.2 
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Figure 8.2 shows the infinitesimal results of 
equation (8.23) and the finite deflection results of equations 
(8.27) and (8.33) for a clamped cylindrical shell for which 
c = 1. Also plotted are Hodge's infinitesimal deflection 
results for a rectangular pressure pulse of the same impulse 
with peak pressure equal to ten times the collapse pressure 
and Jone's infinitesimal and finite deflection results for 
impulsive loading. The approximate results of this section 
are slightly higher than the corresponding results from [8—3]. 
All results are for yield surfaces which circumscribe the 
exact yield surface. 

It should be noted that the analyses of Hodge and Jones 
are valid only for c^ < 6. There is no limitation on c for 
the approximate analysis performed in this section. 

Section 3; Cylindrical Panel Subjected to N on-Uniformly 
Distributed Impulsive Load 

Reference [8-6] reported some experimental results on 
the non- axi symmetric dynamic loading of cylindrical panels. 

Of interest here are some results for aluminum 6061T6 panels 
having clamped longitudinal boundaries and free circumferential 
boundaries which were plastically deformed by detonating sheet 
explosive which covered a portion of the inner surface. The 
midsurface configuration is shown in Figure 8.3. 

Nominal test specimen dimensions are shown in Table 8.1. 
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Table 8.1 



Nominal Cylindrical Panel Dimensions [8-6] 



L 



3.0 inches 



e 



0 



45.0 degrees 



e 



1.5 inches 



e 



33 degrees* 



e 



h 



0.125 inches 



0.091 inches 



0.082 inches 



*Varied slightly with thickness 

The dimensions and 9^ are mode parameters and define 
the regions I, II, and III. 

An approximate, infinitesimal deflection solution has 
been obtained using the general procedure of Chapter 7. 

It was assumed that each region (I, II, III) experienced 
rigid body motion. Region I was assumed to rise vertically 
with a characteristic velocity of w^ at x=L, 9=0.- 

Region II was assumed to rotate rigidly about the support 
line BC (or AD) . 

Region III was assumed to rotate rigidly about the line 
AB (or DC) . 

Because of symmetry about the curves x=L and 9=0, only 
one quadrant (0-x-L; 0-9-9 q) need be analyzed. 
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with the requirement that w and w be continuous at the 
hinge curves between regions, the displacement rates for each 
region are found to be 



Region I 

u = 0 (8.38a) 



V = -Wq sin0 


(8.38b) 


w = Wq cose 




(8.38c) 


Region II 






II 

o 




(8.38d) 


V = n^WoU 


o 

0 
w 

CD 

O 

1 

CD 


(8.38e) 


w = sin ( 00 - 6 ) 


(8.38f) 


, cos 6 

where n = a 

sinOg-ea) 






Region III 

RX w 

u 

^0 


(cos 6 - cos9q) 


(8.38g) 


• X 

V = - — w_ 

=^0 “ 


sine 


(8.38h) 


CD 

tn 

0 

0 

0 

1 

XIX 

11 


(8.38i) 


The positive directions 


of u, V, and w are indicated in 





Figure 8.3. 
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Continuity of w and w also determines the equation of 



the hinge curve C 22 between regions II and III. 



X 



23 ^0 '^a cos8 



sin Oq-B) 



(8.39) 



Equations (8.38 and 8.39) agree with the results of 
Janas [8-7] for the static collapse mode velocities and hinge 
curve of a shallow cylindrical panel before shallow shell 
approximations were made. 

From equations (8.38) the in-surface velocity dis- 
continuities and transverse velocity derivative discontinuities 
are found to be: 

Between regions I and II: 



[u] = 0 



(8.40a) 



[v] = cos(0Q-0^) + sin0^} Wq 

[w'l = 0 

[w*l = {sin0 - n cos (0 -0 )} w 

di 3. U 3. U 



(8.40b) 



(8.40c) 



(8.40d) 



where. 




and 




125 



(8.40e) 



Between regions II and III: 

* R * 

[u] = — (cose - cos6_) w 

Xq u 

[v] = ~ cos(6Q-e) + tan6 sin(6Q-e)} WQ(8.40f) 

w 

[W] = - — cose (8.40g) 

0 

[w*] = - ri^{cos ( e^-e ) - tane sinCe^-e)} Wq (8.40h) 

Between regions I and III: 

* R • 

[u] = — (cose - cose„) w_ 

Xq U U 

[V] = 0 

[W-] = - cose 
^0 

[w*] = 0 

At the clamped boundary; 

[u] = [v] = [w'] = 0 (8.40m) 

[w*] = - (8.40n) 

If the curvature rates for a cylinder are found from 
equations (6.1) and then integrated across a hinge, the 
relative rotation rate components can be shown to be: 



(8.401) 
(8.40j) 
(8.40k) 

(8.401) 
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[W] 



(8.41a) 



and 



Q = 

X 



A _ [v] - [w*] 
“e " R 



( 8 , 



which may be evaluated using equations (8.40). 

With the six-dimensional approximate yield surface 
described in Chapter 7 (Figure 7.2) normality considerations 
for Wq> 0 require the following generalized stress values. 
Between regions I and II: 



M = M ; N = N 
6 0 0 0 



M = N =M_=N. = 0 
X X x6 x6 



( 8 . 



Between regions I and III: 



M = M-; N = N. 
X 0 X 0 






(3 



Between regions II and III; 



M = = 2M - = M- 

X 6 x6 0 



N = N- = 2N - = N- 
X 6 x0 0 



( 8 , 



At the clamped boundary: 



Mg = -Mg 



^x = \ = \0 = \0 = ^0 = ° 



(8 



41b) 



42a) 



42b) 



42c) 



42d) 
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If the appropriate deformation rates (8.40 and 8.41) 
and generalized stress values (8.42) are substituted into 
the internal energy dissipation rate expression (equation 
7.6) and the integration is carried out, it may be shown that 
the dissipation rate for one quadrant of the panel divided 
by the characteristic velocity, w^, and the yield moment. 



M 



0 



IS 



D. 

1 







cos^e (tane -tanO ) } + :^^ {(0 -6 ) 

U U . 0. ^ U 0 



+ cos 6 




sec6 Q+tan0Q 



cos6_ {tan6_-tan6 } 
K 0 U U 0 



+ COS0Q { (tan0Q-tan0^) - j cos0q 



^^^^sec6 +tan0 
0 



sec0Q+tan6 



— ) + tan.9_ sec6 - tan6 sec6 ] } 

U U 0 0 



(8.43) 



Let us call the right hand side of this expression F. 



F . F(L, X^, e^,, e^) 



(8.44a) 



and 



D. 

1 




(8.44b) 
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The external energy dissipation rate (equation 5.1) 



becomes 






Dg = / [(p.-]jw)w - yiiu - yvv] dS 
S quadrant 



(8.45) 



The indicated integration must be performed in each 

region, and the expressions for the pressure and velocities 

within each region depend on the assumed size of the central 

region (Region I) relative to the loaded area (defined by 

I and 0 ) . 
e e 

• • • ^ 

Since u and v are defined in terms of Wq , equation (8.45) 

can be expressed as 



The governing equation of motion is found by equating 
the internal and external energy dissipation rates and has 
the form 



Dg = (P3g - pew„)W|3 



(8.46) 



where 



g - g(L, Xq , 0Q, 0^, 0^) 



and 
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(8.47) 



FMq = ?3g - yewjj 

The upper bound collapse pressure is 




which for a given panel with fixed loading geometry varies 
with the mode shape parameters and 0^. 

The pressure p^ is zero everywhere outside of the loading 
area (Figure 8.3) and g is found by evaluating the integral 



g = / w d S (8.48) 

S ® 

e 



where indicates the loaded portion of the midsurface area. 

Equation (8.48) consists of contributions from each 
region, and the limits of integration within each region 
depend on the parameters Xq and 0^. The form of the equation 
for the values of Xq and 0^ indicated in Figure 8.3, for 
example, would be 



g = 



®a L 
/ / 
0 X, 



R w^ dx d0 -t 



/ 

0 . 



L 

/ 

X 



23 



R Wjj dx d0 



®a L 

if f 

0 X 



0 



/ 

0 . 



^23 

/ } 

x^ 



R dx d0 



(8.49) 
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where the subscript I indicates region I, etcetera. Equation 

(8.39) is used to define X 2 ^ in the limits of integration. 

The optimum values of the mode parameters are those 

which give the lowest static collapse pressure. Because of 

the complicated forms of F and g, parameter optimization was 

accomplished by computing the collapse pressure for discreet 

combinations of and 9 . In all cases it was found that 

the minimum collapse pressure occurred for x.=L and for 0 

U 3l 

somewhat less than 6^. Region I degenerated into a curve 
as shown in Figure 8.4. 

Now e may be determined as 



e = 






6 T 

3. Xj 

[ / / 

0 0 



/ 

0 



0 23 

/ 

0 



] (uii + vv + w\-J) R dx d0 



0 



w^w 



0 L 



/ f (uii + vv + ww)t-_ R dx d9 



0 0 9^ x^3 



II 



(8.50) 



where the mode velocities and accelerations for the first 
set of integrals are those of region III and those of the 
final integral are of region II. Evaluation of this 
expression leads to 



e 

RL 




7 sin20 

a 



2cos0^ sin9 + 0 cos^0») 

0 a a 0 
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Fig, 8.4 
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n R 

+ [7 (0Q-0^) sin0Q + — ^ (sin20Q - sin20^) 



-sin20Q (sin0Q - sin0^) + (Qq”®^) sin0Q cos^0q 



COS0Q 

+ — ^ — (cos20q - cos20^) - 2cos^0q (cos0q - cos0^) 



sec0^ 

- cos^0„ 



n 3 + [0. - - sin(0^-0 ) 

0 sec0_ a 0 a 0 a 



sec0-. 

■" ’'I “"ii^ 



+ sin0. cos^0- (tan0_ - tan0 -0^+0 ) 
0 0 0 a 0 a 



sec0 . 

-2 (0^-0 ) sin0^ + 2cos0- £n 5 — + 2sin20. (sin0rt- 

0 a 0 0 sec0 0 0 



sin0 ) + 2cos20„ (cos0_ - cos0 ) 
a 0 0 a 



- sin20, 



sec0„+tan0. cos 0^ 

Q — r 5 — - — 5 [7(tan^0^ 

sec 0 +tan 0 3 ^ 0 



sec0 



- tan^0 - S,n s^] } 

a sec0 



( 8 . 51 ) 



From the definition of the Martin and Symonds mode factor, 

V, equation ( 8 . 21 ) it can be seen that for constant density, 
y, the denominator of v is peV^. The impulsive velocity, Vq, 
is equal to the explosive impulse divided by the mass of the 
panel portion included within a radial projection of the sheet 
explosive boundaries. 
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"'x' 

' > 



i 








I 



• 'll 

I 



4 




For the degenerate mode shape shown in Figure 8.4, the 



numerator of v is 



yjvj 



®a 

[ / / 

0 L-£ 



®e ^23 

+ f f 

6 L-£ 

a e 



] (yvt 



III 



V. ) Rdx d6 
1 



®e L 

+ / / 

®a ^23 



yV. V. Rdx de 



(8 



where the are the mode velocity components in region A 

and the are the actual impulsive velocity components. 

Performing the indicated integration, one finds that 



J = II + “"la ' < V»a' 

+ COS0 - COS0 - ri (sin^ 0«-co^0 _ ) (0 -0 ) 
e a a"^ 0 Oea 



+ sin0Q 
+ COS0Q 



(cos0 



sec0 



£n 



sec0 



+ 



+ 



COS0 ) + COS0„ 
a 0 

tan0 cos^0 

i + 

tan0 2 

a 



(sin0 - sin0 ) 
a e 

(tan0 -tan0 ) ] } 
e a 



(8 



For the impulsive solution p^ is assumed zero and the 
initial characteristic velocity is assumed to be 



Wq ( t=0) 



vV 



0 



where v 



J 

e 



.52) 



.53) 
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The equation of motion (8.54) becomes 



FM, 



= - 



0 ye 






Motion ceases at t^ when 



WQ(tf) = 0 



and 



yevV, 



t^ = 



f FM, 



™0 t^ 

w = - — ■— + vV t 

0 ye 2 0 



(8.55) 

(8.56) 



(8.57) 



(8.58) 



The final displacement-thickness ratio is found by 
evaluating (8.58) at t^. 



h 2eFMQh 



(8.59) 



Figure 8.5 shows the final central permanent displacement- 
thickness ratios observed during the experiments on the four- 
teen aluminum panels reported in [8-6] plotted as a function 
of a non-dimensional impulse parameter, ip . 






yR V' 



M, 



(8.60) 
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O 0.125 inch 
A 0.091 inch 
D 0.082 inch 



D 




0.618X Circumscribed Yield Surface 



Circumscribed Yield Surface 



Fig. 8.5 
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The same figure contains a plot of points calculated from 
equation (8.59) for panels of the three different nominal 
thicknesses reported in [8.6]. Within plotting accuracy, all 
calculated points fell on the same straight line. The lower 
line of calculated results is based on the approximate yield 
surface which circumscribes the exact six-dimensional Tresca 
yield surface. The upper line of calculated results is based 
on the approximate yield surface whose scale is 0.618 that of 
the circumscribing approximate yield surface. The smaller 
yield surface may not inscribe the exact yield surface, but, 
as was noted in Chapter 7, no approximate yield surface which 
reduces to the two-moment limited interaction yield surface 
and has a scale factor greater than 0.618 can possibly inscribe 
the Tresca yield surface. 

The fact that the smaller yield surface gave better 
results for this problem than the larger yield surface is 
not surprising. A large portion of the hinges has generalized 
stresses corresponding to the corners marked A and B in 
Figure 7.2. Corners of a circumscribed yield surface lie 
outside of the exact yield surface, while the corners of the 
inscribed yield surface may lie on the exact yield surface. 

Section 4: Spherical Cap Subjected to Uniformly Distributed 

Impulsive Load 

As a final example let us consider the case of a deep 
spherical cap with clamped boundary subjected to a transverse 
dynamic pressure load. 
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Giannotti [8-8] reported the results of some experiments 
on clamped spherical caps which were loaded by detonating 
sheet explosive which covered most of the internal surface. 

Of interest here are the results of some experiments on 
6061T6 aluminum caps. Nominal dimensions are listed in 
Table 8.2 and indicated in Figure 8.6. 



Table 8.2 

Nominal Spherical Cap Dimensions [8-8] 





^e 


R 


h 




(degrees) 


(degrees) 


(inches ) 


(inches) 




90 


89 


2.58 


0.115 




90 


89 


2.58 


0.176 




90 


89 


2.58 


0.237 




60 


59 


2.58 


0.117 




45 


44 


2.58 


0.119 




Onat and 


Prager considered 


the static 


collapse problem 


of a deep clamped shell subjected to a uniformly distributed 


pressure [8-9] 


. With an assumed 


velocity 


profile of the 


form 


W^((j),t) = 


Wp(t) (cos4) - COS4 )q) 




(8.61a) 


(1 - COS({)q) 







u = 0 (8.61b) 
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Fig. 8.6 
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where is the transverse velocity at (p=0, and as a conse- 
quence of using the Tresca yield surface, the state of 

generalized stress was found to be membrane (i.e., M. = 

<P 

Mg = 0) except in a narrow region near the boundary, 
cf)* - - 4>g. The angle <p* is defined by the relation 



coscj)* 



2cos(|)g 




( 8 . 62 ) 



derived from the normality requirements of plasticity. 

The upper and lower bound collapse pressures obtained 
by Onat and Prager for shells in the range of and 
listed in Table 8.2 agreed within 15 percent or better. 

The loading used by Giannotti [8-8] appeared sufficiently 
similar to uniform loading to use the Onat and Prager collapse 
profile (equation 8.61). 

For comparison purposes the spherical cap problem was 
solved for infinitesimal deflections and for finite deflections. 
Two infinitesimal solutions with the Onat and Prager velocity 
profile were obtained, one using the Tresca yield surface and 
one using the two-moment limited interaction approximate yield 
surface (Figure 4.1). Lower bounds on the final displacement 
were also obtained using the method of Morales and Nevill 

[8-11] for a perfectly plastic material. 

First let us examine the finite deflection solution. 

In order to simplify the finite deflection analysis, the 

two-direction limited interaction yield surface (Figure 4.2) 
was utilized in place of the exact Tresca yield surface. 
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The strain rate and curvature rate equations (6.1) for 
a spherical cap with the velocity field of equation (8.61) 
reduce to 



w. 



e ^ = 



(J) R(1-cos(})q) 



w sin^cf) 

{cos(|) - coscj). + T T— } 

^0 l-cos(f)- 



(8.63a) 



= 



Wq(cOS(]) - COS(J)q) 
R (1-cos<1)q) 



(8.63b) 



K . = 



o .^2 



Wq C0S(() 



(8. 63c) 



R^ (1-cos({)q) 



The last term of equation (8.63a) is the only finite 
deflection contribution. 

If we assume that the two hinges' rotation rates can be 
represented by average curvature rates, each one averaged over 
one half of the arc RAtf**, in conjunction with a meridional 
strain rate averaged over the arc RAcj)*, it can be shown from 
the normality requirements of plasticity that 



N. „ sin4),, - sincj)* - A(j)* cos({)^ 

(p X _ K U U 

~ h sin<j)Q 

A4>* _ sin2(j)Q - sin2(j)* 



sin<})Q (1 -cos(J)q) 



(8.64a) 



and 
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sin(j)Q - sine})* - A(j)* cos(pQ 
" h sin(f>* 



+ 




^ 4 ,* sin 2 <j>Q - sin 2 (j)* 
~ 4 

sinej)* (l-coscj)^) 



(8.64b) 



where A and B refer to hinges A and B of Figure 8.7. 

The bending moment associated with each hinge may be 
found from the yield surface equations 



_i = 

M 

0 



N' 



±(1 



_ li 



(8.64c) 



n; 



where (-) is used for hinge A and (+) for hinge B when Wq>0. 



For (})-(})* 



% = ^0 



and 



= 0 
4 > 



(8.64d) 

(8.64e) 



Circumferential bending is important only in the same 
region as meridional bending and it can be shown that for 

Wq>0. 



N 

N 



e 

0 




COS(J) ' 



and 



(8.65a) 
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3 








Fig. 8.7 
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4r2 



(1 



cos<p 

°) 2 

COS(J) 



for 4>* - (}) - (f»Q 



while 

= ^0 



(8.65b) 



(8.66a) 



and 

Mg = 0 (8.66b) 

for <})-(})* 



If the velocity profile (equation 8.61), pressure, and 
the generalized stress values (equations 8.65 and 8.66) are 
substituted into the external dissipation rate expression 
for shells of revolution (equation 6.11), and the indicated 
integration performed, it can be shown that 



Rc, Rc- 



P {a, +a_+a .+a^ [ l-{: 



+ P {a„+a 



Rc, 



- 2a., 



RC1C2 



c^ 2 “5 h(l-y) ^*^6 hCl-y) 



Rc, c_ 

2a? ^ } Wq 



+ P {a 



- ^rU 



c ^^"5 hd-y) 6'hd-y) 



- a, wj = 



^8^3 



(8.67) 
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0 







where 



is the Onat and Prager lower bound collapse pressure 

8Mo 



Rh 




Y = 


sin4>* 

sin4>Q 


^1 


sincfjQ - sincj)* - A(J)* coscJjq 
sinct)Q 


^2 


sin2(j)Q - sin2<j)* 

_ 2 4 

sin<j)Q (l-coscj)^) 


^1 


l-cos^** , . . 

= 2 + cos(})q cos(j)* - cos(})q 


^2 


1 ,l-cos»4,* , cos*j,(cos4,*-l) 

R(1-cos4)q) ^3 2 ^ 


^3 


= ^ [i(3cos^(J)Q + cos^(J)*) - 2 cos(J)q coscj)* + cos^cJjq* 


^4 


cos(J)q 

^ " cos^ct)^) - ^ a^] 


^5 


sini^)* (sincj)^ - sine})*) (cos^)* - cos(P^) 
~ 2R (1-cos4>q) 


^6 


= sincj)* (sin(j)Q - sine})*) 
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0 



A 



i 








h sin^<f>Q 
" 8R 

3g = -^ (1 -COS(})q) ^ 

^9 = J (1-cos4>q) ^ 

and 

y is the mass density per unit midsurface area . 

By defining new constants, equation (8.67) can be 
written in the form 

Wq + - P^bgW^ = (8.68) 



where 



®9'"l 


is 


the 


coefficient 


of 






is 


the 


coefficient 


of 


O 

o 


and ^ 9^3 


is 


the 


coefficient 


of 


P w^ on the left hand side of 
c 0 


equation 


(8. 


67) . 










Equation (8.68) has been solved for the impulsive loading 
case by the method of successive approximations [8-10]. 
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The pressure is assumed to be zero and the initial 
characteristic velocity is 



WQ(t=0) = vVq 



(8.69) 



For the partial loading of the spherical cap in 
conjunction with the assumed velocity profile (equation 8.61) 
the Martin and Symonds mode factor becomes 



V = I {l-( 



COS(}) - COS(j), 



1 - COS(}) 



-)M 



0 



(8.70) 



The actual impulsive initial velocity, Vq , is equal to 
the explosive impulse divided by the mass of the spherical 
cap included within a radial projection of the sheet explosive 
boundary. 

Using the second approximation to the solution of 
equation (8.68) one finds 



w^(t) = sinAt + J 2 cosAt + sin^At + cos^At 



+ sin2At + t sinAt + J_ t cosAt + Jq 
5 6 7 o 



(8.71) 



where 



A = /PhZ 
C 2 



^1 = 



-^(1 ^) 
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''2 = e : 



4bjb3 2b3V^Vj 



3b> 






V^v 



■^3 3b, 



0 + 

'=2 



J. = 



3 

3b, 



2v^Vf. b^ 

( + — ) 

b^ 



= 



3b|A 



^6 = 



b^b^A 



J-, = 






and 



b^b 



J„ = - ,— + 



13 



The velocity, w^Ct), can be found by differentiating 
equation (8.71) with respect to time. Motion ceases at the 
time, t^, when w^ = 0. The time t^ is thus the first positive 
root of the transcendental equation obtained by setting w^ 
equal to zero. 

The final central displacement, is obtained by 

evaluating equation (8.71) at time t ^ . 
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Computed results show less than one percent difference 



between the first and second approximations. 

These approximate analytic results are themselves quite 
complicated, due primarily to the presence of finite deflec- 
tion terms. In order to determine whether or not the additional 
complexity is justified, the problem was also solved for 
infinitesimal deflections alone. 

If the same velocity field and the exact Tresca yield 
surface are assumed, and if a hinge is allowed at (})q, the 
dissipation rate expression for shells of revolution 
(equation 6.11) produces the following infinitesimal results 



Equation (8.72) is a linear, ordinary second order 
differential equation. 

For the impulsive case, p^ = 0. Integrating twice with 
respect to time and using the initial conditions that 




(8.72) 



Wq (t=0) = vVq 



and 



Wq ( t=0) = 0 



one finds that 



w^(t) 



P^(a^+2a2+2a^+a^)t^ 




(8.73) 
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Finding from the fact that motion ceases when 
Wo(tf)=0 and substituting into (8.73) leads to the final 
permanent central displacement-thickness ratio 




v‘ 



0 



2P h (a, +2a,+2a .+a_) 
c 1 3 4 7 



(8.74) 



If the Onat and Prager velocity profile is used in 
conjunction with the two-moment limited interaction yield 
surface the following infinitesimal deflection results are 
found 



P 

c 



(l-COS(}>.) ^ 

{ — -H_ + 



h_ 

4R 



(1-cos^4)q) } 




^^) (1-cos(})q) ^ 

(8.75) 



The impulsive final central displacement-thickness ratio 



is 



w 



Of 



y V' 



0 



3P_h(l+^^ 



(8.76) 



2R 



cos 



<}>o' 



Morales and Nevill [8-11] have developed a lower bound 
displacement technique for rigid-plastic structures subjected 
to impulsive loading. The method is based on the use of a 
kinematically admissible velocity field. 

The lower displacement bound obtained from the Onat and 
Prager velocity field and the Tresca yield surface is 
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(8.77) 



w 



Of > 



h 




V 



2 

0 



3P h (a, +2a_+2a ,+a_) 
cl 3 4 7 



The lower bound is seen to be two-thirds the value 
predicted by the Tresca yield surface infinitesimal solution. 

Comparison of experimental and approximate theoretical' 
results are shown in Figures 8.8 through 8.13. These figures 
show final permanent central displacement-thickness plotted 
as a function of impulse parameter X. 



X = 



yR^vJ 

-M^ 



4pR^v2 



(8.78) 



where p is density based on volume. 

Figure 8.13 is a combined plot of the information from 
Figures 8.8, 8.9, and 8.10. This figure shows that 
if one combines the experimental results for the 
individual series of 90 degree caps the complete set of 
results approximates a curve. 

Each figure (except 8.13) has a plot of experimental 
results, Tresca infinitesimal results, two-moment limited 
interaction infinitesimal results, Tresca lower displacement 
bound, and finite deflection results. 

It is immediately obvious that in the range of observed 
displacement-thickness ratios the effect of finite deflections 
on the analytic solution is negligible. 
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Equation (8.73) 

Finite deflection (from 8.71) 
Equation (8.76) 

Experimental results 
Lower bound 

Fig. 8.8 
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Equation ( 8 , 73 ) and finite deflection 
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Experimental results 
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Fig. 8. 9 
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Fig. 8.10 
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Eig. 8.11 
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Equation (8,73) 

Finite deflection (from 8.71) 
Equation (8.76) 

Experimental results 
Lower bound 



Fig. 8.12 
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Fig. 8.13 
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It is also obvious that agreement with experimental 
results is only fair. Particularly noteworthy is the fact 
that some experimental results for the 60 degree caps are 
below the infinitesimal deflection lower bound for a perfectly 
plastic material and all experimental results for the 45 degree 
cap are below the lower bound. 

The 90 degree cap approximate theoretical results can be 
plotted in a manner which show final displacement to be 
essentially independent of thickness for a given initial 
velocity. This behavior is indicative of membrane behavior. 
Indeed bending does have very little influence in the analysis 
of this problem. The experimental results do not show this 
independence of thickness, which indicates that bending may 
be of greater importance than assumed. 

In the interests of simple analysis it was assumed that 
u=0. In order to include greater bending dependence, it may 
be necessary to utilize a more complicated velocity field. 

The simple infinitesimal results obtained for this 
problem may, however, serve as useful engineering estimates 
of shell response. 
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CHAPTER 9 



SUxMMARY AND CONCLUSIONS 

Exact solutions of the response of shells dynamically 
loaded above the static plastic collapse load are quite 
complicated. General approximate theoretical procedures are 
developed herein for arbitrarily shaped shells with arbitrary 
boundary conditions and with arbitrary loading made of rigid- 
perfectly plastic material. The influence of finite trans- 
verse deflections is considered. 

Approximate results were compared with exact analytic 
solutions for complete spheres and circular cylinders and 
also with experimental results for cylindrical panels and 
spherical caps. 

The approximate results agreed exactly with the exact 
analysis for a complete sphere and agreed quite well with 
the exact analysis of a complete cylinder loaded axisymmet- 
rically. The approximate results for partially loaded 
cylindrical panels also agreed favorably with experimental 
results . 

The approximate results for explosive loading of 
spherical caps underpredicted some experimental results and 
overpredicted others. The variation was as great as 60% of 
the predicted results. 

More scatter was observed in the experimental results 
for spherical caps than has been observed in related dynamic 
experiments (e.g. the cylindrical panel results). Some of 
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the spherical cap experimental results are below a lower 
bound on displacement for infinitesimal deflection of a 
perfectly plastic material. 

It appears from the approximate results obtained that 
the influence of finite deflections is quite small if the 
shell geometry requires the development of sizeable membrane 
stresses with infinitesimal deflections. Therefore, it seems 
reasonable to use infinitesimal theory to predict permanent 
displacements for such shells. 

It would be desirable to develop the six-dimensional 
generalized stress space Tresca yield surface at least to 
determine the size ratio between inscribed and circumscribed 
approximate yield surfaces. 

In future experiments concerning the dynamic behavior of 
explosively loaded shells, it would be instructive to deter- 
mine photographically whether or not the assumed impulsive 
initial velocity profile is correct. 

Finally, it would be worthwhile to apply the techniques 
of this thesis to other shell geometries and conditions as 
more experimental results become available. 
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APPENDIX A 



RELATION BETWEEN THE TWO-DIRECTION LIMITED INTERACTION 
YIELD SURFACE AND THE TRESCA YIELD SURFACE 

In Chapter 4 it was stated that the two-moment limited 
interaction yield surface and the two-direction limited 
interaction yield surface (Figures 4.1 and 4.2) circumscribe 
the exact Tresca yield surface and that if each were reduced 
in size to 0.618 times original size they would inscribe the 
Tresca yield surface. 

This appendix contains as a sample calculation the 
limiting condition for the two-direction limited interaction 
yield surface. 

The procedure used to inscribe an approximate yield 
surface within the Tresca yield surface is outlined in 
reference [A-1] . It consists of assuming a state of stress 
corresponding to a vertex of the approximate yield surface 
and then finding the common positive factor by which the 
stresses must be multiplied in order for the vertex to lie on 
the Tresca yield surface. The smallest value of the scale 
factor must be used to inscribe the approximate yield surface 
within the Tresca yield surface. 

Since yield surfaces exhibit symmetry, not every combina- 
tion of vertices has to be examined. Symmetry of the two- 
direction limited interaction yield surface is such (see 
Figure 4.2) that only the first quadrant of each curve need 
be examined. 
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A limiting vertex of the approximate surface is the 
state of stress 




If these values are multiplied by the common scale 
factor ^ and the values n^^ = m^ = 0 and n 2 = ^ are 

substituted into the defining equations of the Tresca yield 
surface [A-1] , the limiting equation is found to be 

m^ - m^ = 1 - (nj^-n2)^ (A.l) 

which becomes 

C = 1 - {A. 2) 

Solving for the positive value of the scale factor yields 

-1 + /s’ 

^ ^ -- = 0.618 (A. 3) 

2 
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